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I.  Introduction  &  Synopsis 


The  response  of  organisms  to  metal  exposure  is  dependent  on  the  chemical  and 
physical  associations  (speciation)  as  well  as  the  concentration  of  the  respective 
metal.  Metal  toxicity  is  regulated  by  the  biogeochemical  environment  in  which 
it  exists,  with  the  metals  physical  form,  kinetic  lability,  and  oxidation  state 
mediating  bioavailability.  With  few  exceptions,  the  uptake  rate  of  a  trace  metal 
by  an  organism  is  largely  dependent  on  its  “free”  or  hydrated  concentration  in 
the  growth  environment.  The  “free”  or  available  metal  is  regulated  by  a  complex 
series  of  competition  reactions  among  aqueous  ligands  and  cell  membrane- 
associated  ligands  that  transport  metals  into  the  organism  (Figure  1). 
Therefore  total  and  operationally  defined  “dissolved”  metal  concentrations  are 
not  necessarily  predictive  of  toxicity  or  potential  risk  to  the  environment. 
Current  regulatory  frameworks  and  metrics  do  not  adequately  address  metal 
speciation  and  addressing  speciation  will  enable  more  realistic  estimations  of 
acute  and  chronic  risk  from  metal  exposures.  The  overarching  goal  of  this 
SERDP  project  was  to  advance  our  understanding  of  metal-ligand  binding  in 
order  to  further  the  development  of  practical  and  predictive  models  of  trace 
metal  bioavailability.  This  goal  is  illustrated  in  Figure  2.  Our  primary 
speciation  focus  was  on  Strong  Metal-Binding  Ligands  and  Colloidal  Phases. 

Figure  1.  Copper  speciation  in  aquatic  environments. 
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Figure  2.  Develop  Practical  Speciation  Procedures,  with  Bioassay 
Grounding  for  the  Prediction  of  Bioavailable  Forms  of  Metals. 


In  this  project  we  isolated  important  metal-ligand  pools  and  characterized  the 
ligands  and  their  metal-binding  properties  in  three  contrasting  marine 
estuaries,  focusing  primarily  on  Copper  (Cu)  and  to  a  lesser  degree,  Zinc  (Zn). 
Chemical  measures  of  metal  speciation  were  linked  to  bioavailability  as 
quantified  in  parallel  laboratory  studies  with  marine  algae.  Multiple 
bioavailability/ toxicity  endpoints  including:  (1)  cellular  budgets  of  trace  metals, 
(2)  molecular  biomarkers,  and  (3)  growth  characteristics,  were  measured. 
Rigorous  trace  metal  “clean”  protocols  (Figure  4)  were  integrated  throughout 
all  components  of  the  study,  including  the  bioassay  studies,  and  state-of-the- 
art  analytical  techniques  were  applied  for  speciation  measurements.  Our 
chemical  speciation  tools  included  electrochemical  methods  (Figure  6), 
principally  voltammetry  (Cathodic  Stripping  Voltammetry  (CSV)  and  Anodic 
Stripping  Voltammetry  (ASV)),  and  kinetic  separations  on  chelating- resins. 
Ultrafiltration  (Figure  5)  at  1  kDa  and  10  kDa  was  also  routinely  applied  to 
fractionate  aquatic  ligand  pools.  For  our  study  systems,  we  chose  three  DoD- 
impacted  marine  estuaries  with  major  contrasts  in  ligand  source,  type,  and 
abundance,  as  well  as  significant  gradients  in  trace  element  levels.  These  are: 
(a)  San  Diego  Bay,  California,  (b)  Norfolk  Harbor,  Virginia,  and  (c)  Cape  Fear, 
North  Carolina.  The  overall  approach  is  diagramed  in  Figure  3. 
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Figure  3.  Project  Approach. 
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Extensive  field  studies  documented  the  existence  of  large  gradients  in  metal 
speciation,  both  between  and  within  systems.  In  San  Diego  Bay,  levels  of 
strong  ligand  [LI]  increase  markedly  from  North  to  South  Bay,  however,  in 
general  concentrations  of  [LI]  and  Cu  are  similar.  In  both  the  Cape  Fear  and 
Norfolk  systems,  levels  of  strong  ligand  are  in  large  excess  of  Cu  levels;  (e.g. 
[LI]  =  lOx  [Cu]),  and  “free”  Cu  levels  are  therefore  extremely  low  (pCu  =  14-15). 
Significant  differences  in  the  “quality”  of  the  organic  matter  are  apparent,  as 
shown  by  large  contrasts  in  normalized  Cu-binding  ligand  levels.  Colloidal 
phases  of  Cu  and  DOM  were  significant,  if  not  the  dominant,  pools  in  each  of 
the  study  systems,  particularly  at  sites  with  high  algal  production  or  where 
there  are  large  inputs  of  terrestrial  carbon.  Cu-binding  ligands  are  also 
predominantly  found  in  the  colloidal  size  fraction,  and  in  both  colloidal  and 
dissolved  phases,  represent  only  a  very  small  fraction  of  DOM  mass.  Kinetic 
separations  on  Chelex  resin  revealed  the  presence  of  large  non-labile  pools  of 
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Cu  in  each  of  the  study  systems,  and  a  close  relationship  was  observed 
between  colloidal  and  non-labile  species.  In  contrast  to  Cu,  a  super-majority  of 
filterable  Zn  in  each  system  was  labile  to  Chelex.  Copper  partitioning  to 
colloids  is  significantly  greater  than  that  to  particles. 

Bioavailability  studies  indicated  that  Cu  uptake  into  cells  can  be  accurately 
predicted  from  levels  of  strong  ligand  [LI].  The  fraction  of  filterable  Cu 
bioavailable  is  greatest  in  the  San  Diego  system.  Bioavailability  is  greatly 
enhanced  by  ultrafiltration,  where  up  to  75%  of  ambient  Cu  is  accessible  after 
ligand  level  reduction.  These  results  strongly  support  efforts  to  develop 
speciation-based  water  quality  criteria.  The  technology  and  models 
underpinning  these  findings  should  be  transferable  to  scientists  within  EPA 
and  DoD. 


Figure  4.  Trace  metal  “clean”  sampling  in  San  Diego  Bay. 
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Figure  5.  Ultrafiltration  of  Cape  Fear  estuarine  water. 


Figure  6.  Cathodic  stripping  voltammetry. 
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II.  Manuscripts  Published  in  Peer-Reviewed  Journals 


The  following  research  manuscripts,  all  published  in  high-impact  factor  peer- 
reviewed  journals,  were  produced  exclusively  from  data  developed  over  the 
course  of  this  study.  SERDP  was  clearly  acknowledged  as  the  funding  agency 
in  each  manuscript.  A  brief  summary  of  our  findings  is  presented  after  each 
manuscript  abstract.  The  full  manuscripts  are  provided  in  the  Appendix  to 
this  report. 


1.  Tang,  D.,  K.  Vang,  D.A.  Karner,  D.A.  Armstrong,  and  M.M.  Shafer.  2003. 
Determination  of  Dissolved  Thiols  Using  Solid  Phase  Extraction  and  HPLC 
Analysis  of  Fluorescently  Derivatized  Thiolic  Compounds.  Journal  of 
Chromatography  A,  998(1-2):31-40. 

Abstract 

A  method  employing  solid-phase  extraction  coupled  with  HPLC  separation  of  thiol- 
monobromobimane  (mBBr)  derivatives  was  developed  and  optimized  to  quantify  dissolved  thiols 
at  concentrations  as  low  as  0.1  n/W  for  glutathione  (GSH)  and  y-glutamylcysteine  (yEC)  in 
natural  waters.  Careful  control  of  the  derivatization  conditions  is  crucial  for  successful  detection 
in  both  the  direct  and  the  solid-phase  extraction  measurements.  The  reducing  reagent,  tri-n- 
butylphosphine  (TBP),  is  needed  for  complete  derivatization.  At  the  optimal  addition  of  TBP 
([TBP]/[mBBr]  =  0.4-1 .6).  Cu,  which  is  present  at  higher  concentrations  than  thiols  in  most 
natural  waters,  does  not  interfere  with  the  determination  of  thiols  under  these  derivatization 
conditions.  Other  soft  metals  such  as  Ag(l)  and  Hg(ll)  are  present  at  levels  typically  less  than  10 
p M  in  natural  waters,  so  their  effects  on  thiol  quantification  should  be  minimal.  The  thiol 
fluorescence  signal  was  totally  suppressed  if  the  mole  ratio  of  TBP  to  mBBr  was  2.6  or  greater. 
Consistent  recovery  of  thiols  standards  in  a  NaCI  solution  (0.5  M)  was  obtained  using  the 
Waters  HLB  reversed-phase  resin,  and  blank  levels  of  GSH  and  yEC  were  extremely  low  (less 
than  0.03  r\M).  The  detection  limits  for  GSH,  yEC  and  phytochelatin-2  (PC-2)  were  0.03,  0.03, 
and  0.06  n/W,  respectively.  Linear  ranges  for  the  thiols  studied  were  excellent,  typically  from  0.1 
to  20  n/W..  This  study  demonstrates  the  capability  of  a  new  method  for  detection  of  low  levels  of 
dissolved  thiols  in  aqueous  systems  using  solid-phase  extraction  coupled  with  HPLC  detection. 
The  method  will  significantly  enhance  the  ability  to  investigate  the  release  of  dissolved  thiols 
from  algae  and  their  role  in  regulating  trace  metal  bioavailability  in  surface. 


2.  Shafer,  M.,  S.  Hoffmann,  J.  Overdier,  and  D.  Armstrong.  2004.  Physical  and 
kinetic  speciation  of  Copper  and  Zinc  in  three  geochemically  contrasting  marine 
estuaries.  Environmental  Science  and  Technology  38(1 4):381 0-381 9. 


Abstract 

The  physical  and  kinetic  speciation  of  Cu  and  Zn  in  three  impacted  marine  estuaries  was 
examined.  Contrasts  in  sources  of  metal-binding  ligands,  solution  chemistry,  and  hydrologic 
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forcing  between  and  within  the  three  study  systems  (Cape  Fear  River  Estuary,  North  Carolina; 
Norfolk-Hampton  Roads-Elizabeth  River,  Virginia;  San  Diego  Bay,  California)  were  exploited  to 
enhance  our  understanding  of  Cu  and  Zn  speciation.  Trace  metal-optimized  tangential  flow 
ultrafiltration  at  1  kDa  nominal  molecular  weight  limit  (NMWL)  was  used  to  fractionate  <0.4  pm 
species  into  colloidal  and  “dissolved”  pools.  Colloidal  species  of  dissolved  organic  matter  (DOM) 
and  copper  were  significant  and  often  the  dominant  pools  in  each  of  the  three  study  systems. 
Characteristic  colloidal  fractions  of  both  DOM  and  Cu  ranged  from  near  70%  of  <0.4  pm 
concentrations  in  Cape  Fear  to  50%  in  San  Diego  Bay.  Colloidal  Cu  and  DOM  were  strongly 
coupled,  and  variability  in  observed  <0.4  pm  Cu  concentrations  was  closely  related  to  the 
concentrations  of  colloidal-associated  metal.  Colloidal  fractions  were  much  smaller  for  Zn  than 
that  of  Cu;  ranging  from  10-30%  in  Cape  Fear  to  less  than  5%  in  San  Diego  Bay,  and  no 
relationship  to  DOM  was  observed.  Kinetic  separations  on  Chelex  resin  revealed  the  presence 
of  large  non-labile  pools  of  Cu  in  each  of  the  study  systems,  with  the  highest  fractions  (70- 
100%)  in  Cape  Fear  and  Norfolk  and  lowest  (30-50%)  in  San  Diego  Bay.  A  close  relationship 
was  observed  between  colloidal  and  non-labile  Cu  species,  implying  slow  reactivity  of  colloidal- 
bound  Cu.  The  fraction  of  filterable  Zn  labile  to  Chelex  averaged  97%,  85%,  and  60%  in  San 
Diego,  Norfolk,  and  Cape  Fear,  respectively.  Anthropogenic  Zn  appeared  almost  exclusively  in 
the  <1  kDa  fraction,  while  anthropogenic  Cu  was  distributed  between  dissolved  and  colloidal 
pools.  Copper  particle-partition  coefficients  (Kd)  followed  the  trend:  San  Diego  »  Norfolk  > 
Cape  Fear  and  were  inversely  correlated  with  DOC  concentrations.  Colloid-based  partition 
coefficients  were  significantly  greater,  in  many  cases  an  order  of  magnitude  greater,  than 
particle-based  partition  coefficients.  The  partitioning  data  suggest  the  presence  of  metal- 
enriched  bacterial-derived  exudates  and/or  discrete  metal  phases  in  colloidal-sized  particles  in 
impacted  regions  of  these  estuaries.  The  strong  relationships  observed  between  Cu  and  DOC 
indicates  that  Cu  partitioning  behavior  over  a  range  of  estuarine  environments  may  be  modeled 
effectively  with  a  limited  set  of  coefficients.  Our  measurements  of  metal  lability  and  size 
distribution  imply  that  the  fraction  of  <0.4  pm  Zn  that  is  likely  to  be  bioavailable  is  greater  than 
that  for  Cu,  especially  in  impacted  regions  of  the  study  systems. 


3.  Tang,  D.,  M.  Shafer,  D.  Karner,  J.  Overdier,  and  D.  Armstrong.  2004.  Factors 
affecting  the  presence  of  dissolved  glutathione  in  estuarine  waters. 
Environmental  Science  and  Technology  38(1 6):4247-4253. 


Abstract 

We  investigated  factors  influencing  the  presence  of  the  thiol  glutathione  (GSH)  in  estuarine 
waters  and  summarize  the  current  state  of  knowledge  about  the  presence  of  GSH  in  surface 
estuarine  waters.  Our  study  addressed  thiol  phase-association,  the  biological  release  of  GSH 
from  a  coastal  diatom  ( Thalassiosira  weissflogii)  in  response  to  short-term  exposures  to  Cu,  and 
the  role  of  copper  (and  other  trace  metals)  in  both  thiol  release  and  preservation  in  surface 
waters.  We  present  new  data  on  GSH  distribution  in  three  diverse  marine  estuaries  in  the 
continental  United  States  (San  Diego  Bay,  Cape  Fear  Estuary,  and  Norfolk  Estuary), 
demonstrating  that  dissolved  GSH  is  present  at  sub-nanomolar  levels,  and  is  preferentially 
partitioned  into  the  ultra-filtrate  fraction  (<1  kDa)  in  comparison  with  dissolved  organic  carbon 
(DOC).  Concentrations  of  GSH  generally  increased  with  increases  in  total  copper  (Cu)  levels, 
although  large  variability  was  observed  among  estuaries.  The  lack  of  a  strong  relationship  of 
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GSH  concentration  across  systems,  sites,  and  season  to  a  single  variable  is  not  surprising  in 
view  of  the  expected  influences  of  biota  species  type,  other  metals  (e.g.,  Mn),  and  concurrent 
and  possibly  opposing  effects  of  Cu  and  DOC  on  GSH  release  and  preservation.  In  30-h 
exposure  experiments,  release  of  dissolved  GSH  from  the  diatom  Thalassiosira  weissflogii  into 
organic  ligand-free  experimental  media  was  a  strong  function  of  added  Cu  concentration.  The 
released  GSH  increased  from  about  0.02  to  0.27  fmol/cell  as  Cu  was  increased  from  the 
background  level  (0.5  nM)  to  310  nM  in  the  modified  Aquil  media.  However,  excretion  of  GSH 
was  lower  (up  to  0.13  fmol/cell)  when  cells  were  grown  in  surface  waters  of  San  Diego  Bay, 
despite  much  higher  total  Cu  concentrations.  Experiments  conducted  in-situ  in  San  Diego  Bay 
water  indicated  that  high  concentrations  of  added  Cu  destabilized  GSH,  while  both  Mn(ll)  and 
natural  colloids  promoted  GSH  stability.  In  contrast,  laboratory  experiments  in  synthetic  media 
indicated  that  moderate  levels  of  added  Cu  enhanced  GSH  stability.  In  summary,  we  have 
shown  that  GSH,  the  major  low  molecular  weight  thiol  in  phytoplankton,  is  released  from 
phytoplankton  in  significant  amounts  and  that  release  is  enhanced  by  Cu.  Given  the  observed 
lability  of  GSH  in  solution  and  our  measurements  documenting  the  presence  of  significant  levels 
of  GSH  in  natural  waters,  it  is  plausible  that  the  dissolved  species  of  GSH  detected  is  the  more 
stable  disulfide  form  (GSSG)  and/or  complexed  with  metals.  The  influences  of  Cu  on  GSH 
stability  are  complex.  Our  experiments  indicated  that  added  Cu  promoted  degradation  of  GSH  in 
a  natural  water  but  enhanced  stability  in  a  synthetic  medium.  These  differences  may  be  due  to 
contrasts  in  GSH  concentration,  GSH/Cu  ratios,  and/or  chemical  composition  of  the  natural 
water  and  synthetic  medium.  Dissolved  organic  matter  may  promote  stability  through 
conjugation  or  enhance  degradation  through  photochemical  reactions. 


4.  Tang,  D.,  M.  Shafer,  D.  Karner,  and  D.  Armstrong.  2005.  Response  of  non¬ 
protein  thiols  to  copper  stress  and  extracellular  release  of  glutathione  in  the 
diatom  Thalassiosira  weissflogii.  Limnology  and  Oceanography  50(2):51 6-525. 


Abstract 

In  this  manuscript  we  present  results  from  bioassay  experiments  using  a  coastal  diatom, 
Thalassiosira  weissflogii,  data  that  document  (1)  the  physiological  response  of  the  diatom  to 
copper  exposure  in  terms  of  photosynthetic  pigments  and  intracellular  thiols,  (2)  the  induction 
and  release  of  essential  thiols  from  the  diatom  cells  and  the  implications  of  these  processes, 
and  (3)  comparability  of  cell  growth  and  excretion  of  metal-complexing  ligands  in  the  metal- 
buffered  and  metal-unbuffered  media.  We  studied  the  dynamic  changes  of  cellular  thiols  and 
the  extracellular  release  of  glutathione  (GSH)  during  growth  of  the  marine  diatom  Thalassiosira 
weissflogii  under  varying  levels  of  copper  (Cu)  addition  in  both  metal  buffered  (with  EDTA)  and 
un-buffered  (without  EDTA)  media.  In  summary,  we  found  that,  with  certain  exceptions,  cell 
growth,  intracellular  thiol  induction,  and  extracellular  GSH  release  were  comparable  at  similar 
Cu  exposure  levels  in  both  the  EDTA-buffered  and  un-buffered  medium.  These  findings 
substantiate  the  premise  that  algae  respond  nearly  exclusively  to  inorganic  species  of  Cu.  In 
both  media,  specific  growth  rates  of  greater  than  1  per  day  were  obtained  at  total  inorganic  Cu 
concentrations  of  less  than  80  nmol  L'1;  however,  at  higher  Cu  levels,  cell  growth  was 
significantly  suppressed.  The  cell  quotas  of  thiols  and  Chi  a  both  decreased  with  growth  time, 
so  that  Chi  a- normalized  cellular  thiol  concentrations  were  more  or  less  conservative,  with 
normalized  values  in  the  range  of  0.5  and  1.5  mmol  GSH  (g-Chl  a)'1.  A  clear  dose-response 
relationship  was  observed  between  PC2:GSH  ratio  and  inorganic  Cu  levels  in  EDTA-un- 
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buffered  medium,  regardless  of  growth  time;  however,  a  more  complicated  pattern  was 
observed  in  EDTA-buffered  medium.  GSH  was  released  from  the  phytoplankton  cells  at 
generally  similar  concentrations  into  both  the  EDTA-buffered  and  un-buffered  growth  media — so 
it  does  not  appear  that  the  presence  of  a  synthetic  metal  chelator  substantially  affects  thiol 
release.  GSH  release  was  closely  related  to  Cu-induced  cell  membrane  damage.  The 
extracellular  GSH  release  rate  was  higher  in  normally  grown  cultures  than  in  growth-limited 
cultures  but  lower  than  in  growth-suppressed  cultures.  There  is,  however,  an  indication  that 
EDTA  possibly  enhanced  the  release  of  GSH,  as  indicated  by  the  higher  deaddive  cell  ratios 
under  EDTA-replete  conditions.  We  also  demonstrated  that  the  release  of  GSH  from  the  diatom 
and  Cu  exposure  are  in  some  sense  decoupled.  For  normally  growing  cells,  about  1.5-2. 2%  of 
cellular  GSH  is  released  in  the  early  exponential  growth  phase  and  9.2%  released  in  the  late 
exponential  phase,  with  the  difference  likely  reflecting  physiological  changes  during  growth 
(mainly  the  decrease  of  cellular  GSH  quota).  At  background  Cu  conditions  (~1  nmol  L'1), 
glutathione  was  excreted  at  an  average  rate  of  0.087  fmol  cell"1  d"1  in  both  media,  but  release 
increased  with  increasing  inorganic  Cu  concentrations.  At  the  elevated  Cu  exposures  (total 
inorganic  [Cu]  100  nmol  L"1),  substantially  greater  amounts  of  GSH  were  released, 
corresponding  to  the  changes  of  the  cell  membrane  integrity.  The  excretion  of  GSH  apparently 
reflects  physiological  conditions  during  algal  growth  rather  than  an  enzymatic  response  of  the 
algae  to  control  trace  metal  speciation  in  the  media.  Thus,  the  Cu-enhanced  release  of 
glutathione  into  ambient  waters  is  probably  an  inadvertent  by-product  of  cell  membrane  damage 
by  Cu  rather  than  a  feedback  mechanism  to  control  the  speciation  (and  toxicity)  of  aqueous  Cu. 
However,  at  background  levels  of  Cu,  glutathione  released  from  algae  could  contribute  a 
significant  portion  of  the  Cu-complexing  ligands  in  oceanic  waters. 


5.  Karner,  D.,  M.  Shafer,  J.  Overdier,  J.  Hemming,  and  W.  Sonzogni.  2006.  An 
algal  probe  for  copper  speciation  in  marine  waters:  Lab  method  development. 
Environmental  Toxicology  and  Chemistry  25(4):  1 106-1 113. 


Abstract 

The  present  experiments  were  designed  to  develop  a  laboratory-based  marine  algal  bioassay  to 
examine  the  bioavailability  of  Cu  with  minimal  use  of  synthetic  ligands,  such  as  EDTA.  Our 
principal  research  objective  was  to  develop  a  robust  bioassay  protocol  for  metal  speciation  in 
marine  waters  using  an  algal  probe.  We  wanted  to  minimize  trace  levels  of  EDTA  in  samples  to 
avoid  compromising  ambient  speciation  but  also  to  allow  chemical  speciation  of  metal-binding 
ligands  excreted  by  the  phytoplankton  cells,  which  may  be  masked  in  the  presence  of  EDTA. 
The  principal  focus  of  our  experiments  was  on  Cu;  however,  some  experiments  with  Zn  also 
were  conducted  (not  reported  here).  We  found  that  Zn  toxicity  was  substantially  lower  than  that 
of  Cu  and  likely  was  irrelevant  in  most  marine  waters.  The  marine  diatom  Thalassiosira 
weissflogii  was  employed  as  our  bioassay  test  organism.  The  first  goal  of  the  present 
experiments  was  to  determine  the  optimal  assay  length,  temperature,  and  light  conditions 
resulting  in  optimal  Cu  concentration  responses.  The  second  goal  was  to  determine  if  algal 
growth  could  be  successfully  maintained  when  EDTA,  Cu,  and  Zn  were  removed  from  the 
experimental  media.  The  final  goal  was  to  establish  if  robust  Cu  concentration  responses  could 
be  achieved  without  EDTA  buffering. 
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Laboratory-based  algal  assays  were  developed  to  explore  the  bioavailability  of  copper  to  the 
marine  alga  Thalassiosira  weissflogii.  The  assay  was  shown  to  be  both  reliable  and 
reproducible  to  increasing  levels  of  Cu.  Growth  consistently  decreased  as  total  Cu  increased  to 
340  nM,  with  a  significant  reduction  in  growth  occurring  at  total  Cu  levels  of  100  nM  or  greater. 
A  calibration  strategy  was  developed  that  avoided  use  of  the  synthetic  ligand 
ethylenediaminetetraacetic  acid  (EDTA)  in  the  Aquil  growth  medium,  thereby  allowing  ambient 
metal  speciation.  In  a  comparison  of  T.  weissflogii  cells  grown  in  Aquil  medium  with  EDTA  to 
medium  containing  no  added  copper,  zinc,  and  less  than  0.003  nM  of  EDTA,  no  significant 
growth  differences  were  observed  after  8  d,  indicating  adequate  stored  nutrients.  A  30-h  assay 
was  selected  as  the  optimal  time  frame  after  examination  of  data  from  concentration-response 
experiments.  Using  65Cu  stable  isotope  additions,  parameters  examined  included  growth, 
chlorophyll  a,  copper  uptake,  phytochelatin  production,  and  dissolved  organic  carbon  excretion. 
The  T.  weissflogii  specific  growth  rates  decreased  from  1.36  d"1  at  pCu  (i.e. ,  the  negative 
logarithmic  concentration  of  free  Cu)  =  8.8  to  0.56  d'1  at  pCu  =  7.8,  whereas  intercellular  copper 
concentrations  increased  from  13.6  to  70.1  fg/cell,  respectively.  Calculated  values  of  the  copper 
concentration  that  caused  a  50%  reduction  in  algal  growth  of  pCu  =  7.7  and  copper  per  algal 
mass  of  625  pg/g  were  established.  Using  an  algal  assay  based  on  EDTA-free  culture  medium, 
along  with  trace-metal  clean  techniques,  the  effect  of  copper  on  T.  weissflogii  and  the  speciation 
of  copper  in  marine  waters  can  be  studied.  The  method  described  has  been  applied  to 
contrasting  marine-estuarine  environments,  and  the  results  from  these  studies  will  be  presented 
in  follow-up  papers. 


The  elimination  of  EDTA  from  growth  media  such  as  Aquil  has  the  potential  to  improve 
phytoplankton-trace  metal  interaction  studies  in  several  important  areas.  Its  removal  not  only 
eliminates  direct  algal  cell  toxicity  associated  with  high  concentrations  of  EDTA  but  also  EDTA- 
enhanced  release  of  glutathione,  an  antioxidant  and  primary  precursor  to  phytochelatin 
production.  The  presence  of  EDTA  severely  compromises  phytoplankton  studies  addressing 
intra-  and  extracellular  complexation  of  metals  by  competing  with  algal-produced  -  excreted 
ligands  (e.g.,  thiols,  phytochelatins)  for  the  target  metals.  The  presence  of  EDTA  in  Aquil  also 
can  bind  essential  metals  required  by  phytoplankton.  As  demonstrated  in  the  present  study,  use 
of  a  large  excess  of  EDTA  in  the  Cu-spike  solution  to  achieve  low  free  Cu  levels  may  have  the 
unintended  consequence  of  binding  other  essential  trace  metals,  such  as  Zn,  and  decreasing 
them  to  a  level  at  which  growth  is  suppressed.  A  recent  study  demonstrated  that  high 
concentrations  of  EDTA  can  change  the  speciation  of  Fe  and  also  raise  organic  carbon  levels; 
therefore,  the  use  of  EDTA  was  not  recommended  when  studying  Fe  limitation  in  phytoplankton. 
For  many  metals  of  concern,  effective  control  of  free  metal-ion  levels  over  a  wide  concentration 
range  is  difficult  to  achieve  with  EDTA,  and  use  of  alternative  synthetic  chelators  can  be  equally 
(or  more)  problematic.  Although  use  of  synthetic  ligands  to  control  free  ion  levels  in  metal- 
bioassay  experiments  has  its  place,  the  present  study  and  other  recent  investigations 
demonstrate  that  the  use  of  EDTA  is  not  necessary  and  may  even  be  problematic.  By  applying 
innovative  trace-metal  clean  methods,  use  of  EDTA  in  trace-metal  bioassay  experiments  can  be 
avoided. 
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6.  Karner,  D.,  M.  Shafer,  J.  Overdier,  J.  Hemming,  and  W.  Sonzogni.  2007.  Roles 
of  dissolved  and  colloidal  DOC  in  regulating  copper  toxicity  to  the  marine  alga 
Thalassiosira  weissflogii.  In  review  at:  Environmental  Toxicology  and  Chemistry. 


Abstract 

The  specific  influence  of  colloidal  (1  kDa  -  0.4  pm)  and  dissolved  (<  1  kDa)  natural  organic 
matter  (NOM)  on  copper  (Cu)  bioavailability  and  toxicity  to  the  marine  alga  Thalassiosira 
weissflogii  was  evaluated  in  water  samples  collected  from  three  estuaries:  Cape  Fear,  NC; 
Norfolk,  VA;  and  San  Diego,  CA  (USA).  Trace  metal  clean  filtration  and  ultrafiltration  techniques 
were  employed  to  separate  samples  into  the  <  0.4  pm  and  <  1  kDa  fractions  (colloid-free)  used 
in  30  h  bioassay  titrations  with  65Cu  additions.  Bioavailability  was  determined  by  measuring  the 
cellular  Cu  quota,  whereas  toxicity  was  indicated  by  a  reduction  in  the  specific  growth  rate 
SGR).  With  increasing  Cu  concentrations,  the  specific  growth  rates  generally  decreased 
(controls  with  no  added  Cu  averaged  1.23  and  1.22  d"1  for  <  0.4  pm  and  <  1  kDa  fractions, 
respectively;  decreasing  to  0.70  and  0.45  d'1,  respectively,  at  highest  Cu  additions)  and  cellular 
Cu  concentrations  generally  increased  (controls  averaged  5  and  4  fg/cell,  respectively,  for  <  0.4 
pm  and  <  1  kDa  fractions,  increasing  to  101  and  129  fg/cell,  respectively,  at  highest  Cu 
additions).  Colloidal  NOM  was  shown  to  exhibit  a  clear  protective  effect  as  average  IC50  values 
from  the  <  1  kDa  fractions  (402  nM  Cu)  were  lower  than  those  from  the  <  0.4  pm  fractions  (492 
nM  Cu).  Partition  coefficients  demonstrate  that  Cu-binding  to  algal  particles  is  greater  at  San 
Diego  relative  to  Norfolk  and  Cape  Fear  and  inversely  related  to  NOM  levels.  When  normalized 
on  an  organic  carbon  basis,  colloidal  organic  matter  in  all  three  systems  was  more  effective  in 
reducing  Cu  toxicity  to  the  algal  cells  than  dissolved  organic  matter.  Samples  collected  nearest 
to  the  bay-mouth  at  all  three  harbors  exhibited  particularly  efficient  colloidal  NOM  sequestration 
of  Cu,  emphasizing  the  importance  of  this  complexation  in  mixing  zones.  We  demonstrate  that 
colloidal  and  dissolved  organic  carbon  can  play  important  roles  in  mitigating  the  bioavailability 
and  toxicity  of  trace  metals  to  marine  phytoplankton  in  estuarine  systems. 

In  summary,  we  have  presented  Cu-titration  data  for  three  marine  estuaries  using  the  marine 
diatom  T.  weissflogii  as  a  probe  of  DOC  characteristics.  We  have  directly  measured  and 
compared  the  influence  of  natural  colloidal  and  dissolved  carbon  on  the  bioavailability  of  Cu  to 
the  alga.  The  calculated  concentrations  that  cause  a  50%  decrease  in  growth  were  consistently 
higher  in  the  <  0.4  pm  fractions  compared  to  the  <  1  kDa  fractions  demonstrating  the  important 
role  colloidal  organic  matter  has  in  buffering  toxicity  of  trace  metals,  such  as  Cu,  to 
phytoplankton  cells.  Copper  partitioning  to  the  diatom  cells,  indicated  by  KD  values,  is  inversely 
related  to  DOC  concentrations  and  follows  the  trend  San  Diego  »  Norfolk  >  Cape  Fear. 
Colloidal  phases  in  both  allochthonous  DOM  dominated  systems  (Cape  Fear)  and  in 
autochthonous  dominated  systems  (San  Diego)  were  better  able  to  protect  alga  from  Cu  toxicity 
than  dissolved  phases  (when  normalized  to  respective  organic  matter  levels).  At  all  three 
harbor  mouth  locations,  colloidal  SGR/DOC  ratios  were  consistently  lower  than  the  dissolved 
SGR/DOC  ratios,  indicating  that  colloidal  phases  are  regulating  Cu  speciation  at  coastal  ocean 
sites  where  low  DOC  levels  are  typically  observed. 
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(Environmental  Chemistry  &  Technology)  -  University  of  Wisconsin-Madison. 
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V.  Field  Campaigns 

Field  studies  were  designed  to  exploit  contrasts  in  sources  and  nature  of  metal 
ligands,  solution  chemistry,  and  hydrologic  forcing  to  enhance  our 
understanding  of  copper  and  zinc  speciation. 

A.  Study  Sites 

I.  San  Diego  Bay 

The  principal  hydrologic  and  geochemical  characteristics  of  San  Diego  Bay 
driving  trace  element  speciation  are:  (1)  Minimal  fluvial  loading,  (2)  Primarily 
autochthonous  dissolved  organic  carbon,  (3)  Generally  high  total  dissolved  Cu 
levels,  and  (4)  A  ligand  gradient  from  North  to  South  Bay. 

Figure  7.  Sampling  Sites  in  San  Diego  Bay. 
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Five  high  intensity  focus  sites  (F)  and  five  non-focus  (NFF)  sampling  sites  were 
established  in  the  San  Diego  Bay  harbor/ estuary  system  (Figures  7  and  8). 
Sampling  at  the  non-focus  sites  was  equivalent  to  that  at  the  focus  sites  except 
that  ultrafiltration  was  not  carried-out.  Coding  and  details  of  site  locations 
are  provided  in  Table  1.  Thirty- three  sites  were  occupied/ sampled  over  the 
course  of  the  study. 

The  San  Diego  studies  were  carried  out  with  close  collaboration  and  integration 
with  Drs  Zirino’s  and  Chadwick’s  group.  The  sampling  strategy  of  the  San 
Diego  Bay  work  was  designed  around  a  residence  time  model  of  the  Bay 
developed  by  Dr  Chadwick.  Sampling  sites  were  positioned  along  a  gradient 
from  short  (North  Bay)  to  long  (South  Bay)  water  residence  time.  In  addition  to 
the  primary  gradient,  sites  which  reflected  end-member  sources  and  highly 
anthropogenically  altered  regions  were  sampled. 

Figure  8.  Aerial  Photo  of  San  Diego  Bay. 
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II.  Norfolk  Harbor 


The  principal  hydrologic  and  geochemical  characteristics  of  Norfolk  Harbor 
driving  trace  element  speciation  are:  (1)  Extensive  fluvial  loading,  (2) 
Autochthonous  +  Terrestrial  dissolved  organic  carbon,  (3)  Large  total  dissolved 
Cu  gradients,  and  (4)  Moderate  DOC  and  ligand  levels. 


Figure  9.  Sampling  Sites  in  Norfolk  Harbor. 
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Figure  10.  Aerial  Photo  of  Norfolk  Harbor. 


Five  high  intensity  focus 
sites  and  five  non-focus 
sampling  sites  were 
established  in  the 
Norfolk/ Hampton 
Roads/ Elizabeth  River 
harbor/estuary  system 
(Figures  9  and  10). 
Sampling  at  the  non-focus 
sites  was  equivalent  to  that 
at  the  focus  sites  except 
that  ultrafiltration  was  not 
carried-out.  Coding  and 
details  of  site  locations  are 
provided  in  Table  1. 
Twenty-two  sites  were 
occupied/ sampled  over  the 
course  of  the  study. 


III.  Cape  Fear  Estuary 

The  principal  hydrologic  and  geochemical 
characteristics  of  the  Cape  Fear  Estuary  driving 
trace  element  speciation  are:  (1)  Fluvial 
dominated,  (2)  Terrestrial  dissolved  organic 
carbon,  (3)  Low  total  dissolved  Cu  concentrations, 
and  (4)  Very  high  DOC  and  ligand  levels. 

Three  high  intensity  focus  sites  and  one  non-focus 
sampling  sites  were  established  in  the  Cape  Fear 
Estuary  system  (Figures  11  and  12).  Sampling  at 
the  non-focus  sites  was  equivalent  to  that  at  the 
focus  sites  except  that  ultrafiltration  was  not 
carried-out.  Coding  and  details  of  site  locations 
are  provided  in  Table  1.  Seven  sites  were 
occupied/ sampled  over  the  course  of  the  study.  A 
sharp  gradient  from  terrestrial  ligands  to 
autochthonous  oceanic  ligands  is  present  from 
sites  CFNC01  to  CFNC03. 


Figure  1 1 .  Aerial  Photo  of  Cape  Fear. 
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Figure  12.  Sampling  Sites  in  the  Cape  Fear  Estuary. 
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B.  Field  Operations  Summary 


Our  primary  field  fractionation  tools,  0.4  pm  filtration,  10  or  1  kilodalton 
ultrafiltration,  chelex  kinetic  separations,  (and  recently  0.02  pm  Anotec 
filtrations)  were  applied  to  all  sites.  Full  scale  ultrafiltration  fractionations  at  a 
molecular  weight  cutoff  of  lkD  were  conducted  at  all  the  focus  sites.  In  San 
Diego  and  Cape  Fear,  splits  of  these  samples  were  distributed  to  both 
Dr.Zirino’s  and  Dr  Skrabal’s  groups  respectively,  for  supplemental 
electrochemical  and  photochemical  experiments.  Chelex  kinetic  separations 
were  performed  at  all  Cape  Fear,  Norfolk,  and,  San  Diego  sampling  sites. 

Bioassay  experiments  utilizing  two  species  of  marine  algae  were  conducted  at 
every  sampling  site  in  the  Cape  Fear,  Norfolk  and  San  Diego  systems.  In 
addition,  at  each  of  the  focus  sites,  detailed  dose-response  studies  were 
carried-out  on  the  bioassay  organisms  by  exposing  the  organisms  to 
progressively  increasing  levels  of  the  65Cu  stable  isotope.  Dose-response 
studies  were  carried-out  on  both  0.4  pm  filtered  samples  and  lkD  ultrafiltered 
samples. 

Samples  for  detailed  electrochemical  speciation  were  obtained,  which  included, 
in  addition  to  the  field  filtered  (0.4  pm)  samples  from  all  sampling  sites,  sub¬ 
samples  from  the  ultrafiltration  fractionations.  Other  samples  collected 
included  those  for  filtered  and  particulate  trace  metals  (including  Cu  and  Zn), 
metal  stable  isotopes,  dissolved  and  particulate  organic  carbon,  sulfide, 
dissolved  and  particulate  thiols,  pigments,  suspended  particulate  matter,  and 
SUVA.  Standard  hydrographic  measurements  were  also  taken. 

In  addition  to  our  prescribed  field  programs,  we  also  implemented  a  series  of 
dedicated  field  experiments  designed  to  follow-up  and  address  emerging  critical 
issues  affected  copper  and  zinc  speciation  in  the  impacted  harbors. 

Table  1  on  the  following  page  summarizes  the  routine  classes  of  samples 
collected  at  focus  and  non-focus  sites  on  a  typical  sampling  campaign. 
Additional  collections  are  made  to  support  dedicated  field  experiments  (e.g. 
incubation  experiments). 
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Table  1.  General  Classification  of  Samples 


Sample  Class 

Focus  Site 

Non-Focus  Site 

Trace  Metals  -  Un filtered  (total) 

X 

X 

Trace  Metals  -  0.4  pm  Filtered 

X 

X 

Trace  Metals  -  Particulate  (>1.0  pm) 

X 

X 

Voltammetry  -  0.4  pm  Filtered 

X 

X 

Voltammetry  -  lkD  Permeate 

X 

— 

Voltammetry  -  lkD  Retentate 

X 

— 

Chelex  -  0.4  pm  Filtered 

X 

X 

Chelex  -  lkD  Permeate 

X 

— 

Bioassay  (Spel)  -  0.4  pm  Filtered 

X 

X 

Bioassay  (Spel)  -  1  kD  Permeate 

X 

Bioassay  (Spel)  -  1  kD  Retentate 

X 

Bioassay  (Spel)  -  Unfiltered 

X 

Bioassay  (Spe2)  -  0.4  pm  Filtered 

X 

X 

Bioassay  (Spe2)  -  1  kD  Penneate 

X 

Bioassay  (Spe2)  -  1  kD  Retentate 

X 

Cu/Zn  Stable  Isotopes  -  0.4  pm  Filtered 

X 

X 

Cu/Zn  Stable  Isotopes  -  Particulate  (>1  pM) 

X 

X 

Organic  Carbon  -  0.4  pm  Filtered 

X 

X 

Organic  Carbon  -  Particulate  (>.7  pM) 

X 

X 

UV  Absorbance  -  0.4  pm  Filtered 

X 

X 

Pigments  by  HPLC  -  Particulate  (>.7  pM) 

X 

X 

Sulfide  -  0.4  pm  Filtered 

X 

X 

Thiol  -  0.4  pm  Filtered 

X 

X 

Thiol  -  Particulate  (>.7  pm) 

X 

X 

SPM  -  Particulate  (>.4  pM) 

X 

X 

Nutrients  -  0.4  pm  Filtered 

X 

X 

Table  2  on  the  following  page  presents  specific  details  on  sampling  site  date 
and  location. 
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Location 

Sampling 

Date  Time 

Site  Code 

Location 

Latitude 

North 

Longitude 

West 

Norfolk 

05-0ct-00 

10:00  AM 

FOCV01 

Upper  Elizabeth  near  Degaussing  Station 

36°  51.834 

076°  19.997 

Norfolk 

05-0ct-00 

2:00  PM 

NRFVA2 

Off  S.  Graney  Island  near  Coast  Guard 

36°  52.663 

076°  20.662 

Norfolk 

05-0ct-00 

4:00  PM 

NRFVA3 

Elizabeth  River  near  Downtown  Norfolk 

36°  51.066 

076°  18.401 

Norfolk 

05-0ct-00 

6:00  PM 

NRFVA4 

Nansemond  River 

36°  50.334 

076°  21.765 

Norfolk 

06-0ct-00 

9:30  AM 

NRFVA5 

Hampton  Roads  near  1-64  Bridge 

36°  58.267 

076°  18.266 

Cape  Fear 

08-0ct-00 

10:30  AM 

CFNC01 

Upper  Cape  Fear  River  at  Wilmington 

34°  10.297 

077°  57.300 

Cape  Fear 

09-0ct-00 

11:00  AM 

CFNC02 

Mid-Lower  Cape  Fear  River 

34°  02.800 

077°  56.509 

Cape  Fear 

10-0ct-00 

10:00  AM 

CFNC03 

Mouth  of  Cape  Fear  River 

33°  55.162 

078°  00.499 

Cape  Fear 

10-0ct-00 

12:30  PM 

CFNC04 

Mid-Upper  Cape  Fear  River 

34°  06.875 

077°  56.087 

San  Diego 

30-Jan-01 

9:00  AM 

SDF01 

Bay  Mouth 

32°  42.930 

117°  13.225 

San  Diego 

31 -Jan-01 

8:15  AM 

SDF02 

Off  Harbor  Island 

32°  43.401 

117°  11.135 

San  Diego 

31 -Jan-01 

1:00  PM 

SDF03 

South  Bay,  Past  Bridge,  Center  of  Channel 

32°  39.861 

117°  08.689 

San  Diego 

01 -Feb-01 

10:00  AM 

SDF04 

Far  South  Bay 

32°  37.396 

117°  07.510 

San  Diego 

02-Feb-01 

9:00  AM 

SDF05 

South  Bay,  Between  Marina  &  Carrier  Dock 

32°  38.883 

117°  08.047 

San  Diego 

29-Jan-01 

3:30  PM 

SDNF01 

Just  South  of  Coronado  Bridge 

32°  40.848 

117°  08.988 

San  Diego 

29-Jan-01 

5:00  PM 

SDNF02 

Off  Shelter  Island  from  Beach  Area 

32°  42.930 

117°  13.224 

San  Diego 

30-Jan-01 

3:30  PM 

SDNF03 

North  of  Coronado  Bridge 

32°  42.159 

117°  10.472 

San  Diego 

31 -Jan-01 

4:30  PM 

SDNF04 

Shelter  Island  Marina 

32°  42.858 

117°  13.861 

San  Diego 

01 -Feb-01 

2:45  PM 

SDNF05 

Off  SPAWAR-Point  Loma 

32°  42.430 

117°  13.945 

San  Diego 

10-M  ay-01 

11:30  AM 

SD0501-A 

Bay  Mouth 

32°  40.923 

117°  13.492 

San  Diego 

10-M  ay-01 

2:00  PM 

SD0501-B 

Just  South  of  Coronado  Bridge 

32°  40.857 

117°  09.089 

San  Diego 

11 -May-01 

9:00  AM 

SD0501-C 

Far  South  Bay 

32°  37.355 

117°  07.401 

San  Diego 

11 -May-01 

11:00  AM 

SD0501-D 

South  Bay,  Between  Marina  &  Carrier  Dock 

32°  38.840 

117°  08.132 

San  Diego 

11 -May-01 

12:45  PM 

SD0501-E 

Off  Harbor  Island 

32°  43.435 

117°  11.086 

San  Diego 

11 -May-01 

2:30  PM 

SD0501-F 

Shelter  Island  Marina 

32°  42.969 

117°  13.768 

Norfolk 

1 7  -  J  u  1  -0 1 

10:00  AM 

NRF0701  -FA 

Hampton  Roads  near  1-64  Bridge 

36°  58.705 

076°  18.555 

Norfolk 

1 7  -  J  u  1  -0 1 

2:30  PM 

NRF0701  -FB 

Off  Graney  Island 

36°  55.801 

076°  21.196 

Norfolk 

1 8-Jul-01 

10:00  AM 

NRF0701-FC 

Lower  Elizabeth  River 

36°  47.448 

076°  18.230 

Norfolk 

1 8-Jul-01 

1:30  PM 

NRF0701-FD 

Upper  Elizabeth  near  Shipyards 

36°  50.489 

076°  17.794 

Norfolk 

1 8-Jul-01 

5:30  PM 

NRF0701  -NFE 

Nansemond  River 

36°  50.380 

076°  21.798 

Norfolk 

19-Jul-01 

10:30  AM 

NRF0701  -NFF 

James  River  near  CR686  Boat  Landing 

37°  08.834 

076°  39.961 

San  Diego 

11 -Sep-01 

10:30  AM 

SD0901-FA 

Bay  Mouth,  Between  1st  &  2nd  markers 

32°  40.383 

117°  13.570 

San  Diego 

12-Sep-01 

3:30  PM 

SD0901-FB 

South  Bay,  Between  Marina  &  Carrier  Dock 

32°  38.883 

117°  08.047 

San  Diego 

1 3-Sep-01 

10:00  AM 

SD0901-FC 

Off  Harbor  Island,  Near  Yellow  Bouy 

32°  43.458 

117°  11.146 

San  Diego 

12-Sep-01 

5:00  PM 

SD0901-NFA 

Just  South  of  Coronado  Bridge 

32°  40.848 

117°  08.988 

San  Diego 

1 3-Sep-01 

12:30  PM 

SD0901  -NFF 

Off  Shelter  Island  from  Beach  Area 

32°  42.603 

117°  13.721 

San  Diego 

25-Feb-02 

4:30  PM 

SD0202-FA 

Bay  Mouth 

32°  40.350 

117°  13.450 

San  Diego 

26-Feb-02 

10:00  AM 

SD0202-FB 

South  Bay,  Past  Bridge,  Center  of  Channel 

32°  39.900 

117°  08.825 

San  Diego 

27-Feb-02 

10:45  AM 

SD0202-FC 

Far  South  Bay 

32°  37.142 

117°  07.281 

San  Diego 

28-Feb-02 

9:30  AM 

SD0202-FD 

Off  Harbor  Island 

32°  43.340 

117°  11.907 

San  Diego 

26-Feb-02 

3:00  PM 

SD0202-NFE 

Shelter  Island,  Mooring  by  Park 

32°  42.603 

117°  13.721 

San  Diego 

27-Feb-02 

3:00  PM 

SD0202-NFF 

Shelter  Island  Marina 

32°  42.965 

117°  13.746 

San  Diego 

28-Feb-02 

1:00  PM 

SD0202-NFG 

North  of  Coronado  Bridge 

32°  41.996 

117°  10.074 

Cape  Fear 

1 6-Apr-02 

1:30  PM 

CF-0402-FA 

Mouth  of  Cape  Fear  River 

33°  56.390 

077°  58.804 

Cape  Fear 

1 6-Apr-02 

9:30  AM 

CF-0402-FB 

mid-Lower  River 

33°  59.133 

077°  56.902 

Cape  Fear 

1 6-Apr-02 

3:30  PM 

CF-0402-C 

mid-Upper  River 

34°  06.875 

077°  56.087 

San  Diego 

14-M  ay-02 

8:30  AM 

SD0502-FA 

Bay  Mouth 

32°  40.383 

117°  13.570 

San  Diego 

14-M  ay-02 

10:30  AM 

SD0502-FB 

Off  Harbor  Island 

32°  43.458 

117°  11.146 

San  Diego 

13-May-02 

2:00  PM 

SD0502-FC 

Far  South  Bay 

32°  37.170 

117°  07.204 

San  Diego 

13-May-02 

3:30  PM 

SD0502-FD 

South  Bay,  Past  Bridge,  Center  of  Channel 

32°  39.677 

117°  08.640 

San  Diego 

14-M  ay-02 

11:45  AM 

SD0502-FE 

Shelter  Island  Marina 

32°  42.967 

117°  13.765 

Norfolk 

1 3-Oct-02 

4:00  PM 

NRF1002-FA 

Hampton  Roads  near  1-64  Bridge 

36°  58.638 

076°  18.450 

Norfolk 

1 3-Oct-02 

7:00  PM 

NRF1002-FB 

James  River  near  CR686  Boat  Landing 

37°  04.974 

076°  39.435 

Norfolk 

14-Oct-02 

11:00  AM 

NRF1002-FC 

Lower  Elizabeth  River 

36°  47.482 

076°  18.190 

Norfolk 

14-Oct-02 

5:45  PM 

NRF1002-FD 

Upper  Elizabeth  near  Shipyards 

36°  50.383 

076°  17.732 

Norfolk 

1 3-Oct-02 

4:45  PM 

NRF1002-NFE 

Off  Graney  Island 

36°  55.859 

076°  21.549 

Norfolk 

18-Jun-03 

9:30  AM 

NRF-0603-FA 

Hampton  Roads  Outside  1-64  Bridge 

36°  58.717 

076°  17.704 

Norfolk 

17-Jun-03 

9:45  AM 

NRF-0603-FB 

Southern  Elizabeth  River  near  1-64  Bridge 

36°  45.135 

076°  17.577 

Norfolk 

17-Jun-03 

1:45  PM 

NRF-0603-FC 

Upper  Elizabeth  near  Shipyards 

36°  50.485 

076°  17.809 

Norfolk 

16-Jun-03 

5:30  PM 

NRF-0603-SFD 

James  River 

37°  00.184 

076°  32.290 

Norfolk 

18-Jun-03 

11:30  AM 

NRF-0603-NFG 

Off  Graney  Island 

36°  55.970 

076°  22.285 

Norfolk 

17-Jun-03 

12:00  PM 

NRF-0603-SFH 

Lower  Elizabeth  River 

36°  47.491 

076°  18.191 

Page  22  of  79 


22 


VI.  Speciation  Studies 

A.  Summary 


Colloidal-species  of  organic  carbon  and  copper  are  significant,  if  not  the 
dominant  phase,  in  each  of  the  three  study  systems  (Table  3).  Colloidal-phase 
zinc  levels  are  much  lower.  Organic  carbon  and  Cu  concentrations  are  highly 
correlated;  however  Zn  concentrations  show  no  relationship  to  DOC.  Higher 
colloid  levels  and  fractions  are  typically  observed  at  sites  with  high 
autochthonous  carbon  production  or  where  there  are  large  inputs  of  terrestrial 
carbon.  Cu-binding  ligands  are  found  predominantly  in  the  colloidal  size 
fraction,  and  in  both  colloidal  and  dissolved  phases,  represent  only  a  very  small 
fraction  of  DOM  mass.  Total  Zn-binding  ligand  levels  are  significantly  lower 
than  Cu-ligand  levels  and  colloidal  phases  are  of  only  minor  importance. 


Table  3:  Speciation  Summary 


[DOC]  pM 
%  Colloidal 

DOC 

DOC 

200-1000 

40-75 

200-500 

30-60 

65-170 

20-45 

[Cu]  nM 

Cu 

6-20  (12) 

12-60  (25) 

4-120  (35) 

[Zn]  nM 

Zn 

4-25  (15) 

6-250 

4-300  (110) 

%  Colloidal 

Cu 

45-80 

40-50 

45-50 

%  Colloidal 

Zn 

15-40 

6-12 

<5 

%  Chelex  NL 

Cu 

70-100 

70-100 

40-80 

%  Chelex  NL 

Zn 

18-30 

0-15 

1-3 

[Li] 

Cu 

70-120 

60-200  (80) 

5-90  (40) 

[Li] 

Zn 

5-40 

2-25 

0-5 

Bioavailability 

Cu 

Low 

Moderate 

High 

Bioavailability 

Zn 

Moderate 

High 

High 

Chelex  NL  =  Metal  not  retained  on  chelating  column  (a  measure  of  reactivity). 


Predicted  bioavailability  of  Cu  follows  the  trend:  San  Diego  >  Norfolk  >  Cape 
Fear.  Bioavailability  of  Zn  is  predicted  to  be  high  in  both  San  Diego  and  Norfolk 
and  lower  at  Cape  Fear. 
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B.  Cape  Fear 


Colloidal  species  of  organic  carbon  dominate  this  tidal  estuary  [Figure  13]. 
Nearly  three-quarters  of  the  total  “dissolved”  organic  carbon  (DOC)  is  colloidal 
in  size  at  the  terrestrial  end-member  site,  dropping  to  65%  at  the  ocean  end- 
member  sampling  site.  The  phase  association  and  size  fractionation  of  copper 
is  highly  correlated  with  that  of  (DOC);  in  fact  sites  CFNC01  and  CFNC03  show 
nearly  identical  Cu  fractionation  to  that  of  DOC.  Riverine  input  of  terrestrial 
carbon  drives  the  DOC  dynamics  (and  Cu-complexation)  in  the  upper  reaches 
of  the  estuary. 


Figure  13. 
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Figure  14.  Distribution  of  Cu  and  Zn  between  dissolved  (<1  kD),  colloidal 
(1  kD  -  0.4  jim),  and  particulate  (>0.4  jim)  phases  at  our  study  sites  in  the 
Cape  Fear  Estuary  [*  =  ultrafiltration  not  performed;  filterable  =  <0.4  pm]. 
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A  very  high  degree  of  consistency  is  observed  in  metal  levels  and  phase 
partitioning  at  a  given  site  [Figure  14].  Concentrations  of  both  metals  drop 
significantly  as  one  moves  downstream  from  Wilmington  to  the  mouth  of  the 
Cape  Fear  River.  This  decline  is  greater  for  Zn  (70  nM  to  15  nM)  than  for  Cu 
(24  nM  to  10  nM).  While  colloidal  forms  of  Cu  are  the  dominant  phase, 
colloidal  species  of  Zn  are  much  less  important,  with  only  minor  contributions 
at  the  river  mouth.  In  contrast,  particulate  forms  of  Zn  are  the  dominant 
phase,  with  particulate  forms  of  Cu  contributing,  but  much  less  so.  The 
dissolved  (<  1  kDa)  fraction  of  Cu  is  typically  less  than  30%  of  total  metal  at  all 
sites.  Trace  metal  and  DOC  Levels  in  Cape  Fear  can  be  modeled  with  a  two 
end-member  mixing  model  (Cape  Fear  River  &  Coastal  Seawater)  [Figure  15]. 
DOC  and  metal  in  <1  kDa  phases  exhibit  only  minor  non-conservative 
behavior. 

Figure  15.  Two  end-member  mixing  model  for  Cape  Fear. 
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Cu-complexing  ligand  at  Cape  Fear  [Figure  16]  is  found  to  be  predominantly 
associated  with  colloidal  sized  phases  at  both  upper  and  mid-estuary  sites.  An 
input  of  <1  kD  sized  ligand  in  the  ocean  end-member  site  results  in  a  nearly 
50:50  distribution  between  colloidal  and  dissolved  phases. 
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CFNC03 

Total  Cu  Ligand  =  113  nM 

Dissolved  =  8% 
Colloidal  =  92% 


CFNC02 

Total  Cu  Ligand  =  68.3  nM 

Dissolved  =  40% 
Colloidal  =  60% 


CFNC01 

Total  Cu  Ligand  =  98.1  nM 


Dissolved  =  49% 
Colloidal  =  51% 


Figure  16.  Phase  distribution  of  Cu-complexing  ligand  at  three  sites  at 
Cape  Fear. 
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A  typical  speciation  modeling  result  for  the  Cape  Fear  system  is  presented  in 
Figure  17.  Strong  ligand  is  in  great  excess  of  total  Cu,  all  (99.99%)  Cu  is 
bound  to  the  strongest  ligand  (LI),  and  the  free  [Cu2+]  is  very  low  (pCu  =  14.5), 
or  0.00005%  of  total  copper. 


Figure  17.  Typical  speciation  modeling  output  for  Cu-complexing  ligand 
at  Cape  Fear. 
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C.  Norfolk  -  Hampton  Roads 


The  Norfolk  system  is  more  complex  both  hydrographically  and  geochemically 
than  the  Cape  Fear  system,  however,  colloidal  phases  remain  critically 
important  [Figure  18].  Several  geochemically  contrasting  riverine  sources  of 
Cu-complexing  ligand  are  evident,  including  the  James  River  (with  a  large 
proportion  of  <  1  kD  DOC  and  Cu)  and  the  Elizabeth  River  (with  a  higher 
proportion  of  colloidal  DOC  and  Cu).  As  in  the  Cape  Fear  system,  colloidal  Cu 
levels  are  strongly  related  to  colloidal  carbon  levels;  however,  carbon  levels  are 
significantly  lower,  and  the  nature  of  carbon  more  variable,  so  the  relationship 
is  not  as  tight.  Site  NRF-A  shows  evidence  of  coagulation/ aggregation  of  DOC 
and  Cu-ligand  associated  with  aggregate  precursors. 

Figure  18.  Distribution  of  Cu  and  DOC  between  dissolved  (<1  kD)  and 
colloidal  phases  (1  kD  -  0.4  pm)  at  selected  focus  sites  at  Norfolk. 
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Figure  19.  Distribution  of  Cu  and  Zn  between  dissolved  (<1  kD),  colloidal 
(1  kD  -  0.4  pm),  and  particulate  (>0.4  pm)  phases  at  our  study  sites  in  the 
Norfolk  Estuary  [*  =  ultrafiltration  not  performed;  filterable  =  <0.4  pm]. 
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A  very  high  degree  of  consistency  is  observed  in  metal  levels  and  phase 
partitioning  at  most  sites  [Figure  19]  (e.g.  compare  four  years  of  sampling  at 
Hampton  Rds,  and  four  years  of  sampling  at  the  Upper  Elizabeth  sites). 
Concentrations  of  both  metals  at  the  ocean  end-member  site  (Hampton  Rds) 
are  comparable  to  those  observed  at  the  ocean  end-member  site  at  Cape  Fear, 
however,  river  and  harbor  concentrations  of  metals  are  significantly  greater  in 
the  Norfolk  system  than  at  Cape  Fear  (40-80  nM  versus  20-25  nM  Cu  and  150- 
400  nM  versus  60-70  nM  Zn.  Colloidal  forms  of  Cu  are  a  co-dominant  phase  at 
most  sites,  however,  colloidal  species  of  Zn  are  much  less  important,  with  only 
minor  contributions  at  all  sites  studied.  In  contrast,  particulate  forms  of  Zn  are 
the  dominant  phase  at  all  sites  impacted  by  the  James  River,  including  the 
ocean  end-member.  Particulate  forms  of  Zn  at  the  Elizabeth  River  sites  are 
significant  but  not  dominant.  Significant  (20-40%)  particulate  contribution  to 
total  Cu  levels  is  also  observed.  Dissolved  (<  1  kD)  Cu  levels  are  typically  less 
than  50%  of  total  Cu  levels  at  all  study  sites,  and  a  similar  dissolved  fraction  is 
measured  at  James  River  impacted  sites  for  Zn.  However  in  the  Elizabeth  River 
system  dissolved  fractions  of  Zn  are  much  greater  (>80%). 

Colloidal  sized  (1  kD  -  0.4  pm)  Cu- 
binding  ligands  typically  represent 
50-60%  of  total  ligand  in  the 
Norfolk  system  (Figure  20).  Very 
high  levels  of  ligand  are  measured 
at  site  NRF-C,  the  Elizabeth  river 
end-member  (this  river  drains  a 
large  wetland  and  in  doing  so 
acquires  high  levels  of  DOC). 


Figure  20.  Phase  distribution 
of  Cu-complexing  ligand  at  two 
sites  in  the  Norfolk  system. 
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A  typical  speciation  modeling  result  for  a  Norfolk  site  with  relatively  high  levels 
of  both  total  Cu  and  DOC  is  presented  in  Figure  21.  The  general  copper 
speciation  pattern  is  intermediate  between  that  determined  at  Cape  Fear  and 
that  in  San  Diego  Bay.  Total  Cu  levels  are  approaching  strong  (LI)  ligand 
concentrations,  and  therefore  a  small  fraction  of  Cu  is  bound  to  weaker  L2 
ligands.  However,  little  or  no  Cu  is  associated  with  ASV  ligands.  Free  Cu2+  is 
still  very  small  however  (pCu  =  13.2),  though  an  order  of  magnitude  higher 
than  at  Cape  Fear  when  expressed  as  a  percentage  of  total  Cu.  Modeling 
predicts  a  54%/46%  distribution  between  colloidal  and  dissolved  phases, 
nearly  identical  to  the  ratio  actually  measured  by  ultrafiltration,  lending  strong 
support  to  that  approach. 

Figure  2 1 .  Typical  speciation  modeling  output  for  Cu-complexing  ligand 
at  Norfolk. 
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D.  San  Diego  Bay 


At  the  other  end  of  the  gradient  from  high  DOC  fluvial/ terrestrial  dominated 
systems  to  low  DOC/authothonous  dominated  systems  is  San  Diego  Bay. 
Fluvial  inputs  are  minimal  and  DOC  levels  are  nearly  an  order-of-magnitude 
lower.  DOC  and  ligand  levels  are  regulated  by  a  balance  between  production 
of  authothonous  carbon  and  loss  through  sedimentation/ coagulation  and 
hydrologic  flushing.  Colloidal  carbon  fractions  increase  from  20-25%  at  Bay 
mouth  to  40-50%  in  South  Bay,  a  trend  consistent  with  increasing  water 
residence  time  (Figure  22).  Colloidal  Cu  levels  are  high  (40-55%)  and  are  not 
simply  related  to  bulk 
colloidal  carbon,  suggesting 
the  presence  of  a  very 
strong  ligand,  but  at  low 
concentrations.  The 

observed  trends  in  colloidal 
Cu  fraction  also  reflect  the 
parallel  increases  of  ligand 
and  total  Cu  as  one  moves 
from  Bay  mouth  to  South 
Bay. 


Figure  22.  Distribution 
of  Cu  and  DOC  between 
dissolved  (<1  kD)  and 
colloidal  phases  (1  kD  - 
0.4  jim)  along  a  transect 
from  South  Bay  (SDF04) 
to  Bay  Mouth  (SDFOl)  in 
San  Diego  Bay. 
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Figure  23.  Distribution  of  Cu  and  Zn  between  dissolved  (<1  kD),  colloidal 
(1  kD  -  0.4  pm),  and  particulate  (>0.4  pm)  phases  at  our  study  sites  in  San 
Diego  Bay  [*  =  ultrafiltration  not  performed;  filterable  =  <0.4  pm]. 
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As  observed  in  our  two  other  study  systems,  metal  levels  and  phase 
partitioning  at  most  sites  are  quite  consistent  on  a  seasonal  and  annual  basis 
[Figure  23]  (e.g.  compare  four  samples  each  at  Far  S.  Bay,  Mid  S.  Bay,  and 
Harbor  Is.  sites  -  (the  0501  data  are  an  exception  however)).  Concentrations 
of  Cu  at  the  ocean  end-member  site  (Bay  Mouth)  are  more  variable  (ranging 
from  5-20  nM)  and  typically  lower  than  at  Norfolk  and  Cape  Fear.  This 
reflects  both  the  large  contrast  and  sharp  gradient  in  metal  concentrations 
between  harbor  and  ocean  and  varying  tidal  mixing  conditions.  Zinc 
concentrations,  like  those  of  Cu,  are  also  variable  at  the  ocean  end-member 
site;  however  levels  are  comparable  to  those  measured  at  Norfolk  and  Cape 
Fear.  Copper  and  Zn  concentrations  in  San  Diego  Bay  proper  are  similar  to 
those  measured  in  the  Norfolk  system  (40-60  nM  Cu;  150-200  nM  Zn),  but 
again  significantly  higher  than  at  Cape  Fear  (20-25  nM  Cu  and  60-70  nM  Zn. 
Colloidal  forms  of  Cu  are  a  co-dominant  phase  at  most  sites,  however,  colloidal 
species  of  Zn  are  much  less  important,  with  only  minor  contributions  at  all 
sites  studied.  Unlike  in  the  two  other  study  systems,  particulate  forms  of  Zn 
are  only  a  relatively  minor  component  of  the  total  metal.  The  one  exception  to 
this  “rule”  are  samples  collected  in  May  of  2001,  where  very  large  particulate 
fractions  of  Zn  (and  Cu)  were  measured  -  coinciding  with  a  major  algal  bloom 
in  the  Bay.  Particulate  contributions  to  total  Cu  levels  were  more  significant 
(15-30%)  than  that  measured  for  Zn.  Dissolved  (<  1  kD)  Cu  levels  are  typically 
less  than  40%  of  total  Cu  levels  at  all  study  sites,  however  dissolved  fractions 
of  Zn  are  typically  much  greater  (>85%). 


Except  for  the  ocean  end-member  site  (where  <  1  kD  Cu-binding  ligands 
dominate),  colloidal  phases  completely  dominate  the  Cu-binding  ligand  size 
distribution  (Figure  24).  Between  70  and  80%  of  total  strong-ligand  binding 
capacity  in  San  Diego  Bay  proper  is  associated  with  colloidal  species.  This  is 
in  contrast  to  most  sites  in  the  other  study  systems,  where  strong-ligand 
binding  capacity  more  closely  tracked  DOC  levels.  Note  that  strong  ligand 
levels  double  from  the  Bay  Mouth  (18  nM)  to  North  Bay  (36  nM),  and  double 
again  from  North  Bay  to  South  Bay  (73  nM) 
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Figure  24.  Phase  distribution  of  Cu-complexing  ligand  at  five  sites  in 
San  Diego  Bay,  representing  a  transect  from  Bay  Mouth  (SDF01)  to  South 
Bay  (SDF04). 
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A  typical  speciation  modeling  result  for  San  Diego  Bay  Mouth  Site  is  presented 
in  Figure  25.  Total  Cu  levels  are  similar  to  those  at  Cape  Fear,  however  the 
ligand  levels  are  much  lower,  and  therefore,  the  strongest  ligands  are 
saturated.  Some  Cu  is  bound  to  weaker  ligands  (L2),  a  small  fraction  even  to 
very  weak  (ASV)  ligands.  Free  [Cu2+]  is  much  higher  than  at  Cape  Fear  (pCu  = 
12.5). 


Figure  25.  Typical  speciation  modeling  output  for  Cu-complexing  ligand 
in  San  Diego  Bay  -  Bay  Mouth. 
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A  typical  speciation  modeling  result  for  San  Diego  Bay  -  South  Bay  Site  is 
presented  in  Figure  26.  Levels  of  strong  ligand  and  Cu  are  both  greater  than 
at  the  Bay  mouth.  A  large  majority  of  the  Cu  is  bound  to  the  strongest  ligands 
(likely  algal-sourced),  and  much  less  is  bound  to  weaker  ligands  than  at  the 
Bay  mouth. 

Figure  26.  Typical  speciation  modeling  output  for  Cu-complexing  ligand 
in  San  Diego  Bay  -  South  Bay. 
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E.  Filterable  Copper  -  Strong  Ligand  Comparison 

Total  levels  of  strong  copper-binding  ligand  are  in  large  excess  of  Cu  levels  in 
both  the  Cape  Fear  and  Norfolk  systems  (Figure  27).  At  Cape  Fear,  strong 
ligand  levels  range  from  80  to  over  270  nM  (modal  value  of  -100  nM)  and 
filterable  Cu  concentrations  are  <  20  nM.  At  Norfolk,  strong  ligand  levels  range 
from  60  to  nearly  400  nM,  with  filterable  Cu  concentrations  in  the  range  of  20- 
50  nM.  However  in  San  Diego,  concentrations  of  ligand  and  copper  are  similar, 
and  many  instances,  Cu  levels  exceed  strong  ligand  (Figure  27). 
Bioavailability  of  Cu  is  predicted  to  be  significantly  greater  under  the 
conditions  measured  in  San  Diego  in  comparison  to  other  study  systems. 

Figure  27.  Strong  copper  ligand  concentrations  compared  with 
filterable  copper  levels  at  Cape  Fear,  San  Diego  Bay,  and  Norfolk. 


Page  39  of  79 


39 


F.  Chelex  Lability  Summary 

Uptake  of  copper  and  zinc  onto  a  chelating  resin  (Chelex-100)  was  used  as  a  measure 
of  the  lability  of  metal  in  the  aqueous  samples.  Fundamentally  this  is  a  metric  of  the 
dissociation  kinetics  (lability)  of  metal-ligand  associations  and  therefore  complements 
the  voltammetry  studies  which  are  equilibrium  based.  It  must  be  stressed  that  kinetic 
lability  is  not  necessarily  correlated  with  the  stability  constant  (equilibrium-derived)  -  it’s 
possible  for  a  very  strong  metal-ligand  association  to  exhibit  fast  dissociation  kinetics. 
We  specially  prepared  Chelex  columns  and  operated  them  under  conditions  that  would 
be  expected  to  retain  only  the  free  and  very  labile  metal  ion.  This  work  is  detailed  in  our 
publication:  Shafer,  M.,  S.  Hoffmann,  J.  Overdier,  and  D.  Armstrong.  2004. 

Physical  and  kinetic  speciation  of  Copper  and  Zinc  in  three  geochemically 
contrasting  marine  estuaries.  Environmental  Science  and  Technology 
38(1 4)  .3810-381 9.,  and  only  a  brief  mention  of  this  sizable  effort  is  presented  here. 

Figure  28.  Chelex  non-labile  Cu  versus  colloidal  Cu  in  the  3  systems. 
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A  very  strong  relationship  exists  between  Chelex  non-labile  Cu  (i.e.  the  Cu  not 
retained  on  the  columns)  and  colloidal  copper  (Figure  28),  indicating  that  Cu 
associated  with  colloids  is  likely  not  labile  and  therefore  less  bioavailable.  The 
lines  defining  the  relationships  for  San  Diego  and  Cape  Fear  pass  through  the 
origin,  suggesting  that  all  non-labile  copper  is  represented  by  colloidal  phases. 
However,  at  Norfolk,  a  significant  y-intercept  is  measured  (~10  nM  -  Figure  28) 
signifying  the  presence  of  non-labile  Cu  species  in  the  <  1  kD  dissolved 
fraction.  The  presence  of  strong  synthetic  ligands  such  as  EDTA  or  NTA, 
sourced  from  waste  water  treatment  plant  discharges,  in  the  Norfolk  samples, 
could  explain  this  finding. 

A  comparison  of  the  measured  physical  and  kinetic  speciation  of  Cu  and  Zn  at 
Bay  Mouth  sites  at  Cape  Fear  and  San  Diego  is  presented  in  Figure  29. 
Filterable  (dissolved  +  colloidal)  Cu  levels  are  50%  higher  (10  versus  6.5  nM)  at 
San  Diego;  however,  the  colloidal  fraction  is  significantly  greater  at  Cape  Fear. 
The  labile/ non-labile  distribution  of  Cu  mirrors  that  of  the  dissolved-colloidal 
fractionation. 


Figure  29.  Comparison  of  physical  and  kinetic  speciation  of  Cu  and  Zn 
at  Bay  Mouth  Sites  at  Cape  Fear  and  San  Diego. 
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Filterable  Zn  concentrations  are  comparable  in  the  two  systems,  as  is  the 
distribution  between  dissolved  and  colloidal  phases  -  both  systems  exhibiting  a 
much  greater  Zn  dissolved  fraction  than  that  for  Cu.  The  greater  dissolved 
fraction  for  Zn  is  mirrored  in  the  Zn  lability  data,  with  lability  fractions 
approaching  90%  in  San  Diego.  Though  labile  Zn  fractions  at  Cape  Fear  are 
significantly  greater  than  those  for  Cu,  a  substantial  non-labile  fraction  is  still 
measurable. 


VII.  Ligand  &  Dissolved  Organic  Matter  Character  (DOM) 

Significant  differences  in  the  “quality”  of  the  organic  matter  are  apparent,  as 
shown  by  contrasts  in  carbon  normalized  Cu-binding  ligand  levels  (Figure  30) 
and  Specific  UV-Absorbance  -  SUVA-254  nm  (Figure  31). 

Normalized  dissolved  strong  ligand  levels  in  San  Diego  range  over  a  factor  of  2 
(from  110-210  pM/M).  In  Norfolk,  the  overall  dissolved  levels  are  higher,  but  a 
factor  of  2  range  in  normalized  ligand  abundance  is  also  observed  [from  200- 
400  pM/M],  A  much  larger  range  is  measured  in  the  Cape  Fear  system  (60- 
450  pM/M)  resulting  primarily  from  dilution  of  strong  ligand  by  very  high  DOC 
at  the  river  end-member  sites.  Colloidal  phases  in  all  three  systems, 
particularly  at  sites  with  high  algal  production,  are  dramatically  enriched  in 
Cu-binding  ligand  in  comparison  with  dissolved  phases.  This  colloidal 
enrichment  is  most  evident  in  San  Diego  Bay,  the  system  with  the  greatest 
autochthonous  carbon  production. 

It’s  clear  from  these  normalized  data,  and  very  important  to  mention,  that  the 
strong  Cu-binding  ligand  at  all  sites  and  systems  represents  only  a  very  small 
fraction  of  the  DOM  (irregardless  of  specific  binding  site  stoiciommetry).  At 
typical  values  of  100-200  pM/M  strong  ligand  abundance  corresponds  to  only 
100-200  PPM  in  the  carbon.  Even  in  the  highly  enriched  colloidal  fractions 
strong  ligand  abundance  rarely  exceeds  1000  PPM  (0.1%). 


SUVA-254  data  demonstrate  the  predominantly  terrestrial  character  of  the 
DOM  in  the  Norfolk  and  Cape  Fear  systems  (higher  SUVA  values  are  indicators 
of  increasing  aromatic  content  of  the  DOC).  At  Cape  Fear,  the  ultrafiltration 
retentates  (higher  molecular  weight  fractions)  have  a  greater  aromatic  content 
compared  to  the  permeates  (1  and  10  kDa).  However  in  San  Diego,  the 
opposite  is  true,  again  indicating  a  remarkably  different  pool  of  DOM  than  that 
found  in  the  Cape  Fear  system.  SUVA  values  increase  from  Bay  Mouth  to 
South  Bay  in  San  Diego,  likely  reflecting  the  longer  hydraulic  residence  time 
and  more  extensive  DOM  processing  in  South  Bay. 
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Figure  30.  Comparison  of  organic  carbon  normalized  strong  Cu  ligand 
levels  between  dissolved  (1  kDa  permeate)  and  colloidal  (1  kDa  retentate) 
molecular  weight  fractions  in  our  three  study  systems. 
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Figure  31.  Comparison  of  Specific  UV  Absorbance  (SUVA-  254  nm) 
between  small  (1  kDa  permeate)  and  large  (1  kDa  retentate)  molecular 
weight  fractions  in  our  three  study  systems. 
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Algal  pigment  levels  show  a  strong  correlation  with  Cu-binding  ligand 
concentrations  (Figure  32),  suggesting  that  strong  ligands  are  derived  from 
algal  exudates.  In  San  Diego  the  data  (slopes)  suggests  that  strong  ligand  is 
solely  derived  from  algal  production,  whereas  on  the  East  Coast  a  baseline 
contribution  from  non-algal  sources  is  also  indicated. 

Figure  32.  Strong  Cu-binding  ligand  concentration  -  relationship  to 
algal  pigments  (Fucoxanthin  -  a  diatom  accessory  pigment)  at  study  sites. 
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VIII.  Copper  and  Zinc  Partitioning  to  Particles  and  Colloids 


Our  unique  sampling  and  physical  fractionation  approach  gave  us  the  ability  to 
directly  examine  metal  partitioning  to  colloidal-sized  suspended  material.  The 
findings  of  this  facet  of  the  larger  project  have  been  published  in  a  peer- 
reviewed  manuscript  [Shafer,  M.,  S.  Hoffmann,  J.  Overdier,  and  D.  Armstrong. 
2004.  Physical  and  kinetic  speciation  of  Copper  and  Zinc  in  three  geochemically 
contrasting  marine  estuaries.  Environmental  Science  and  Technology 
38(1 4)  :3810-3819\,  and  therefore  only  a  brief  summary  is  presented  here. 

Summary:  Nearly  without  exception  the  colloidal  phases  exhibit  an 
enhanced  affinity  for  copper  (refer  to  Figures  33  and  34).  Colloids  are  enriched 
in  Cu  (on  a  pg  g-1  basis  relative  to  particles),  and  show  significantly  greater 
partition  coefficients.  Zinc  partitioning  to  colloids,  however,  is  generally  less 
than  or  comparable  to  that  of  particles  (Figure  35).  Dramatic  and  reproducible 
trends  are  apparent  between  study  systems.  For  example  where  Cape  Fear  Cu 
particle  partition  coefficients  were  low  and  relatively  invariant  (4,000-6,000), 
Cu  colloid  partition  coefficients  are  high,  with  clear  trends  in  moving  from  the 
upper  estuary  (101,000  L/Kg)  to  the  lower  estuary  (340,000  L/Kg).  Levels  (pg 
g1)  of  Cu  in  both  colloids  and  particles  are  nearly  a  factor  of  10  higher  in 
Norfolk  than  at  Cape  Fear,  and  levels  at  San  Diego  are  an  order-of-magnitude 
greater  again  than  at  Norfolk.  In  San  Diego  Bay,  colloid  partition  coefficients 
routinely  exceed  one-million.  The  unique  nature  of  the  suspended  phases, 
coupled  with  low  DOC  levels,  results  in  extremely  high  colloid-partition 
coefficients  in  San  Diego  Bay.  These  findings  have  clear  and  dramatic 
implications  for  modeling  the  fate  of  Cu  in  these  three  impacted-systems. 

The  quantity  of  Cu  per  unit  mass  of  suspended  particulate  matter  can  be  used 
as  an  indicator  of  ecosystem  quality.  Our  data  (Figure  33)  shows  dramatic 
and  regular  differences  in  the  magnitude  of  this  descriptor  between  our  three 
primary  study  systems.  At  Cape  Fear  levels  (pg  g1)  are  very  low,  ranging  from 
7  in  the  upper  estuary  to  2  in  the  lower  estuary.  In  the  Norfolk  system  levels 
are  nearly  an  order-of-magnitude  higher  ranging  from  30  to  70.  A  remarkable 
degree  of  consistency  is  observed  between  our  October  2000  and  July  2001 
Norfolk  data  sets.  In  San  Diego  Bay,  Cu  content  of  particles  increases  nearly 
an  order-of-magnitude  above  that  determined  in  Norfolk  (or  lOOx  greater  than 
Cape  Fear).  Though  the  Bay  mouth  site  (ocean  end  member)  was  low  (10-30, 
depending  upon  season),  it  was  not  unusual  for  South  Bay  sites  to  exhibit 
levels  of  >200-400.  One  site  within  Shelter  Island  exhibited  an  extraordinary 
value  of  6300  (0.63%  Cu)  -  which  is  suggestive  of  the  presence  of  discrete 
particles  of  Cu  within  the  water  column. 
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Figure  33.  Copper  content  of  particles  and  colloids  in  the  three  study 
systems. 


Page  47  of  79 


47 


Figure  34.  Copper:  Particle  and  colloid  Partition  Coefficient. 
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Particle-partition  coefficients,  though  not  strictly  thermodynamic  constants, 
are  frequently  used  to  express  the  net  “affinity”  of  a  particle  for  a  metal.  We 
have  determined  these  Kd  values  for  all  our  sites  and  not  unexpectedly  the 
trends  observed  follow  those  determined  for  the  micrograms  per  gram 
descriptor  (above).  Cape  Fear  Kd’s  are  very  low  (4,000-6,000  L/Kg),  reflecting 
both  low  Cu  levels  on  particles  and  competition  from  high  DOC  levels.  Norfolk 
Kd  levels  are  typically  in  the  range  of  20,000  to  40,000,  a  magnitude  not 
atypical  for  many  moderately-impacted  systems.  However,  in  San  Diego  Bay 
KdS  were  nearly  always  greater  than  50,000,  with  several  sites  in  the  hundreds 
of  thousands.  It’s  apparent  that  the  nature  of  the  suspended  particles  in  San 
Diego  Bay  coupled  with  low  DOC  levels,  result  in  a  strong  preference  of  Cu  for 
particulate  phases. 

Cu  Kd-colloid  is  strongly  predicted  in  all  three  study  systems  from  DOC  levels 
(Figure  36).  The  corresponding  DOM-Cu  Kd-particle  relationship  is  significant 
but  not  as  strong  as  the  colloid  relationship. 

Figure  36.  Relationships  between  DOC  and  Colloid  and  Particle 
Partition  Coefficients  for  Copper. 
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IX.  Thiol  Ligand  Production  and  Excretion 


We  carried-out  extensive  and  pioneering  studies  of  thiol  ligand  production  and 
loss  in  estuarine  waters  as  a  major  focus  effort  in  the  overall  project.  This 
effort  was  in  recognition  of  the  fact  that  thiols  likely  play  a  significant  role  in 
sequestering  (binding-ligand)  trace  elements  (particularly  Cu)  in  aquatic 
systems,  and  that  aquatic  organisms  (especially  phytoplankton)  are  major 
players  in  the  generation  and  cycling  of  thiols.  The  cartoon  (Figure  37)  below 
illustrates  processes  affecting  thiol  ligand  production  and  excretion. 
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Figure  37.  Thiol  ligand  production,  excretion,  and  reactions. 

We  quantified  pools  of  speciated  thiols  in  natural  algal  populations  (suspended 
solids)  and  in  the  estuarine  waters.  We  also  directly  measured  thiol  excretion 
by  algae  in  laboratory  experiments  run  in  parallel  with  the  Cu  and  Zn  bioassay 
studies.  Large  scale  field  incubation  studies  were  conducted  in  San  Diego  Bay 
to  ascertain  thiol  excretion  rates  of  indigenous  algae.  Dedicated  thiol 
excretion/ stability  experiments  were  carried-out  on  several  occasions  over  the 
study  period. 


These  studies  produced  three  publications  in  peer-reviewed  journals: 
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Tang,  D.,  K.  Vang,  D.A.  Karner,  D.A.  Armstrong,  and  M.M.  Shafer.  2003. 
Determination  of  Dissolved  Thiols  Using  Solid  Phase  Extraction  and  HPLC 
Analysis  of  Fluorescently  Derivatized  Thiolic  Compounds.  Journal  of 
Chromatography  A ,  998(1 -2): 31 -40. 

Tang,  D.,  M.  Shafer,  D.  Karner,  J.  Overdier,  and  D.  Armstrong.  2004.  Factors 
affecting  the  presence  of  dissolved  glutathione  in  estuarine  waters. 
Environmental  Science  and  Technology  38(1 6) -.4247-4253. 

Tang,  D.,  M.  Shafer,  D.  Karner,  and  D.  Armstrong.  2005.  Response  of  non¬ 
protein  thiols  to  copper  stress  and  extracellular  release  of  glutathione  in  the 
diatom  Thalassiosira  weissflogii.  Limnology  and  Oceanography  50(2) : 516- 
525. 

The  reader  is  directed  to  the  above  manuscripts  (and  synopses  in  the  front  of 
this  report)  for  details  of  these  studies.  Just  a  couple  findings  are  graphically 
presented  below: 


PC-2  in  Thai,  NC0204 


Figure  38.  Phytochelatin  (PC-2)  induction  by  Cu  in  Thallasiosira  (Aquil 
growth  media  (Aq)  ad  Cape  Fear  waters. 

Page  52  of  79 


52 


A  strong  linear  relationship  is  observed  between  Cu  burden  in  the  algae  (as 
well  as  aqueous  Cu  exposure  concentrations)  and  induced  thiol  (phytochelatin 
-  PC-2)  levels  in  the  algae  (Figure  38).  Measured  thiol  induction  slopes  were 
similar  in  defined  culture  media  and  field  site  waters  (Cape  Fear  sites  A  and  B). 


Figure  39.  Dissolved  glutathione  in  San  Diego  Bay 


Field  data  (Figure  39)  confirm  the  predictions  of  laboratory  studies.  Increasing 
levels  of  dissolved  thiol  (glutathione  -  GSH)  are  found  to  be  correlated  with 
increasing  Cu  levels.  The  higher  levels  in  May  compared  with  the  February 
data  set  is  consistent  with  greater  primary  production  observed  during  the  May 
sampling  period. 
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Figure  40.  Glutathione  (GSH)  Release  from  Thallasiosira  in  Aquil  Media 
(GSH)  and  San  Diego  Waters  (FX). 


Increasing  amounts  of  thiol  (glutathione)  are  observed  to  be  excreted  into  the 
growth  media  upon  exposure  to  increasing  levels  of  Cu  (Figure  40).  Less  thiol 
is  excreted  in  the  field  samples  compared  with  the  Aquil  calibration  samples 
which  illustrate  the  protective  effect  of  natural  colloids  on  thiol  induction. 
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X.  Bioavailability  Studies 


Summary:  We  developed  robust  trace  metal  bioassays  using  two  species  of 
marine  phytoplankton  (Thalassiosira  weissflogi,  and  Emiliania  huxleyi ).  In  an 
effort  to  provide  sub-lethal  toxicity  markers,  and  to  support  biotic  ligand 
modeling,  we  examined  a  wide  range  of  endpoints,  including  intra-  and  extra¬ 
cellular  metal  uptake  (levels /burden  in  cells),  phytochelatin  induction  (metal 
stress  proteins),  as  well  as  growth  inhibition.  The  bioassay  methods  developed 
provide  key  values  for  the  levels  of  critical  ligands  (on  a  site  specific  basis) 
protective  of  organism  toxicity  -  i.e.  the  metal  buffering  capacity  of  the  system 
calibrated  to  the  test  organism.  This  is  accomplished  by  titrating  (in  discrete 
increments)  the  ligands  with  added  metal  using  the  test  organism  as  the 
biosensor.  We  also  independently  determined  the  concentrations  of  various 
classes  of  ligands  with  electrochemical  (voltammetric)  techniques  and 
compared  the  biological  endpoints  to  chemical  ones.  This  approach  supported 
one  of  our  principal  project  goals:  i.e.  to  couple  chemical  measures  of  ligand 
quantity  and  strength  to  various  biological  endpoints.  The  bioassay  protocols 
that  we  developed  in  this  project  are  described  in  a  peer-reviewed  publication: 
Karner,  D.,  M.  Shafer,  J.  Overdier,  J.  Hemming,  and  W.  Sonzogni.  2006.  An  algal 
probe  for  copper  speciation  in  marine  waters:  Lab  method  development. 
Environmental  Toxicology  and  Chemistry  25(4) :1 106-1 11 3. 

A.  Brief  Outline  of  Bioassay  Protocols 

i.  Facility 

Algae  are  cultured,  and  most  experiments  performed,  in  the  bio-monitoring 
facility  of  the  University  of  Wisconsin  -  Madison  State  Laboratory  of  Hygiene. 
Experiments  are  carried  out  in  specially  constructed  HEPA-filtered, 
polycarbonate  enclosures,  placed  within  constant  temperature,  walk-in 
incubators.  Sample  processing  is  restricted  to  several  class- 100  HEPA  clean 
benches,  and  trace  metal  clean  technique  is  applied  throughout  all  steps. 

ii.  Algal  Species 

Two  common  marine  phytoplankters  were  evaluated:  (1)  Thalassiosira 
weissflogi,  a  centrales  diatom.  Cells  are  10-12  pm  in  dimension,  and  (2) 
Emiliania  huxleyi,  a  prymnesiophyte  with  CaCCU  coccoliths.  Cells  are  4-8  pm 
in  dimension. 
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Thalassiosira  weissflogi 


Centrales  Diatom  with  Si02 
frustules.  10-12  [am. 
Widely  Distributed 


Emiliania  huxleyi 

Prymnesiophyte  with 
CaC03  coccoliths.  4-8  pm. 
Very  Widely  Distributed 


iii.  Calibration 

Marine  algae  are  cultured  in  an  extremely  well-defined  synthetic  seawater 
media  (modified  -  Aquil).  Aquil  media  is  “chelexed”  to  remove  traces  of 
contaminant  metals  and  filter  sterilized.  (The  Chelex-100,  a  chelating  ion- 
exchange  resin,  is  subject  to  an  exhaustive  10  day  clean-up  protocol  to  remove 
trace  of  metals  and  “loose”  ligand).  Once  sterilized,  the  media  is  handled  under 
axenic  conditions.  Initially,  the  algae  are  grown  in  full  Aquil  media,  which 
incorporates  trace  levels  of  many  micro  nutrients,  including  Cu  and  Zn.  Algae 
are  harvested  in  late  log  phase  from  the  complete  media  and  transferred  to 
Aquil  media  devoid  of  added  Cu,  Fe,  Zn,  and  EDTA.  Timing  is  such  that  after 
6-7  days  in  deficient  media,  growth  dilution  provides  a  stable  and  low 
background  level  of  Cu  and  Zn  in  the  cells,  which  can  then  be  harvested  for 
actual  experimentation.  The  period  from  initial  spiking  of  cells  into  full  Aquil  to 
harvesting  from  deficient  Aquil  in  preparation  for  an  actual  exposure  typically 
spans  two  to  three  weeks. 


Calibration  curves  are  prepared  in  deficient  Aquil  media,  by  the  addition  of  nM 
levels  of  target  metal.  The  metal  is  added  in  the  form  of  a  stable  isotope  (65Cu, 
68Zn).  Prior  to  preparation  of  the  actual  exposure  flask,  the  acid  washed 
polycarbonate  vessels  are  pre-equilibrated  with  an  aliquot  of  Aquil  spiked  with 
metal  to  the  specific  designated  exposure  level.  The  pre-exposure  Aquil  is 
discarded  after  24  hours  of  equilibration.  Dependent  upon  the  algal  species, 
levels  of  spiked  total  metal  range  from  0  to  200  nM,  with  6-7  concentration 
points.  Chemical  modeling  allows  us  to  accurately  determine  the  speciation  of 
Cu  or  Zn  in  the  Aquil  calibration  flasks,  including  levels  of  free  ion. 
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Experiments  are  initiated  by  spiking  into  the  exposure  flasks  the  specially- 
prepared  algal  cells,  typically  about  one  hour  after  the  metal  additions.  For 
Thalassiosira  we  target  for  an  initial  cell  concentration  of  4,000  cells  per  ml;  for 
Emiliania,  15,000  cells  per  ml.  For  a  300  mL  Aquil  exposure  volume,  this  cell 
concentration  is  provided  by  an  addition  of  3-7  mL  of  the  mid-log  phase 
deficient  stock.  Our  standard  exposure  period  is  30  hours  under  an  18/6  light 
cycle.  Additions  are  staggered,  accounting  for  post-exposure  processing  time, 
so  that  variation  in  actual  exposure  periods  among  various  flasks  is  minimized. 

Figure  41.  Thallasiosira  Cu  calibration  curve. 


A  typical  calibration  curve  in  deficient  Aquil  is  shown  above  in  Figure  41. 
Cells  were  inoculated  to  an  initial  density  of  4000  cells  per  mL,  with  an 
exposure  period  of  30  hours.  The  points  and  associated  regression  model  line 
relate  free  Cu  levels  in  the  growth  media  (nM)  to  intracellular  Cu 
concentrations  (fg  per  cell).  Excellent  dose  response  linearity  is  indicated.  Also 
shown  (vertical  bars)  is  the  macro  response  (cell  number  as  determined  by  flow 
cytommetry)  of  the  organism  to  increasing  free  Cu  levels.  A  toxic  “threshold”  of 
approximately  7-8  nM  free  Cu  is  indicated,  with  a  growth  EC50  of  13  nM 
calculated. 
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iv.  Post-Exposure  Processing 

Post-exposure  sub-sampling  is  conducted  under  HEPA  clean  benches,  and 
incorporates  a  series  of  carefully  executed  and  earmarked  filtrations  adapted  to 
the  specific  analyte.  Each  and  every  exposure  flask  undergoes  these  tests  (and 
there  could  by  hundreds  in  a  typical  field  sampling  event).  These  are 
summarized  below: 

1 .  Intra-  and  Extra-  cellular  Trace  Metal  content  of  algal  cells 

-  The  algal  suspension  is  filtered  onto  a  25  mm  pre-cleaned  Teflon 
(1-2  pm)  membrane,  held  in  an  all  Teflon  vacuum  filtration 
apparatus. 

-  To  differentiate  intra-  and  extra-  cellular  associated  metal,  the 
harvested  cells  are  rinsed  with  a  buffered  solution  of  EDTA, 
followed  by  rinses  with  ultra-clean  synthetic  seawater.  Both  the 
filter  and  the  rinses  are  frozen  (-20  °C). 

-  The  filtrate  (after  UV-oxidation  and  buffering)  is  analyzed  for 
metals  either  by  solid-phase-chelation  inductively-coupled  plasma 
mass  spectrometry  (SPC-ICP-MS)  or  high  resolution  sector-field 
ICP-MS  (HR-ICP-MS). 

-  The  filter-borne  algal  cells  are  solubilized  using  an  automated 
microwave-assisted  acid  digestion  in  sealed  Teflon  bombs. 

-  The  algal  digestate  is  analyzed  for  metals  by  solution  ICP-MS. 

2.  Thiol  (including  Phytochelatin)  content  of  algal  cells 

-  The  algal  suspension  is  filtered  onto  a  25  mm  glass-fiber  filter 
(GFF-0.7  pm),  held  in  an  all  Teflon  filtration  apparatus. 

-  The  filter  is  placed  into  a  small  LDPE  vial,  shock  frozen  in  liquid 
nitrogen,  and  then  transferred  to  a  -80  °C  freezer. 

-  Thiols  are  extracted  at  60  °C  using  methanesulfonic  acid  in  a 
borate  buffer,  followed  by  homogenizing  (Teflon  pestle)  and 
centrifugation  to  prepare  a  representative,  clear,  extract. 

-  Thiols  in  the  extract  are  derivitized  with  mBBR  (a  fluorescent  tag) 
in  the  presence  of  TBP  (a  disulfide  reducing  agent). 

-  The  derivitized  thiols  are  quantified  by  HPLC  using  a  dual  solvent 
gradient  of  0.08%  TFA  in  MQ  and  100%  acetonitrile. 

-  Thiols  quantified  include:  Phytochelatin  2  (PC2),  Phytochelatin  3 
(PC3),  glutathione,  cysteine,  glutamylcysteine,  and  acetylcysteine. 

3.  Pigment  content  of  algal  cells 

-  The  algal  suspension  is  filtered  onto  a  25  mm  glass-fiber  filter 
(GFF-0.7pm),  held  in  an  all  Teflon  filtration  apparatus. 

-  The  filter  is  folded,  sealed  in  foil,  placed  into  a  -20  °C  freezer. 

-  The  pigments  are  extracted  from  the  cells  into  a  90% 
acetone / water  matrix. 

-  Pigments  are  quantified  by  HPLC,  using  both  fluorescence  and 
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UV-VIS  (diode  array)  detection. 

-  Typically  the  following  compounds  are  detected  and  quantified: 
fucoxanthin,  diadinoxanthin,  lutein,  (3-carotene,  chlorophyll  a, 
chlorophyll  b,  alloxanthin,  zeaxanthin,  pheophytin,  and 
pheophorbide. 

4.  Cellular  Mass  concentrations 

-  The  algal  suspension  is  filtered  onto  a  tared,  25  mm 
polycarbonate  membrane  (1.0  pm),  held  in  an  all  Teflon  filtration 
apparatus. 

-  The  filter  is  folded,  placed  into  a  40  °C  oven,  and  then  dried  to 
constant  weight.  Filters  are  re-weighed  on  a  micro-balance  to  a 
significance  level  of  0.01  mg. 

5.  Cell  Number  concentrations 

-  An  aliquot  of  the  algal  suspension  is  gently  centrifuged  to  achieve 
a  cell  pre-concentration  factor  of  5. 

-  Cell  are  enumerated  (and  sized)  on  a  Coulter  EPICS  XL  flow- 
cytometer.  Three  separate  aliquots  of  the  concentrated  suspension 
are  measured. 

Other  measurements  that  are  performed  on  a  less  frequent  basis  are  (a)  trace 
metal  content  of  the  filtrate  to  complete  a  mass  balance  (b)  thiol  content  of 
filtrate  to  look  for  excretion  products  (c)  Cathodic  Stripping  Voltammetry  (CSV) 
and  Anodic  Stripping  Voltammetry  (ASV)  to  assess  ligand  release  by  algae. 

The  protocols  outlined  above  allow  us  to  relate  the  free  metal  level  in  the 
growth  media  to  cellular  pools  of  various  constituents,  including  the  trace 
metals. 

v.  Field  Sample  Analysis 

Field  samples  are  handled  in  a  manner  nearly  identical  to  that  described  above 
for  the  calibration  curves,  with  a  few  key  differences: 

a.  Obviously  the  field  sample  (0.4  pm  filtrate,  1  and  10  kD  permeate, 
ultrafiltration  concentrate)  is  substituted  for  the  Aquil  media.  Algae  spiking 
levels,  however,  remain  similar,  with  target  concentrations  of  4,000  and  15,000 
cells /ml  depending  upon  species. 

b.  Metal  stable  isotope  spiking  levels  are  distributed  over  a  significantly  larger 
concentration  range  to  account  for  the  protective  effect  of  natural  aquatic 
ligands  (e.g.  Cu  spikes  at  certain  Norfolk  sites  range  up  to  500  nM). 

Complexation  capacity  and  ECso  are  determined  from  the  inflections  and  slope, 
respectively,  of  the  dose-response  curve.  The  bioavailable  metal  is  probed  from 
the  duplicate  or  triplicate  samples  receiving  no  metal  addition.  The  metal 
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content  of  these  cells  represents  that  which  the  algae  can  acquire  from  the 
natural  environment  -  i.e.  bioavailable.  The  cellular  pool  is  equated  to  the 
labile  pool  in  the  sample  using  the  calibration  relationships  established  with 
each  batch  of  exposures.  The  labile  solution  pool  of  metal  can  be  expressed  as 
a  fraction  of  total  measured  Cu  in  the  sample  to  facilitate  comparison  among 
samples. 


B.  Selected  Outcomes  from  Bioassay  Experiments 

Figure  42.  Comparison  of  Cu  uptake  from  0.4  pm  filtered  samples  by 
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Thallasiosira  at  three  study  sites. 

A  comparison  of  Cu  uptake  from  0.4  pm  filtered  samples  by  Thallasiosira  at 
three  study  sites  (Norfolk  Site  A,  July  2001  -  at  entrance  to  Hampton  Roads; 
Norfolk  Site  C,  July  2001  -  Elizabeth  River  end-member,  San  Diego  Site  A, 
September  2001,  North  Bay,  at  entrance  to  Harbor)  is  shown  in  Figure  42. 
Copper  intra-cellular  uptake  is  measured  in  fg  per  cell,  after  30  hour 
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exposures.  Dose-response  is  typically  linear  at  each  site,  and  uptake  patterns 
are  consistent  with  relative  ligand  levels  between  sites  (e.g.  Norfolk-C  ligand 
levels  are  2.5x  greater  than  at  Norfolk-A  which  results  in  a  dramatic  decrease 
in  uptake  due  to  competition  from  “dissolved”  ligands). 


Figure  43.  Copper  uptake  (fg/cell  Cu)  and  Cell  Growth  (cells/ ml)  of 
Thallasiosira  from  1  kD  permeate  at  Norfolk  Site  D. 
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Dose-response  (Growth  and  Cu-uptake)  of  Thallasiosira  to  Cu  in  1  kD 
ultrafiltration  permeate  from  Norfolk  Site  D  (similar  to  Site  C  in  Figure  42 
above)  is  shown  in  Figure  43.  The  protective  effect  of  the  colloidal  ligands  has 
been  removed  and  toxicity  is  evident  at  much  lower  Cu  concentrations.  The 
cellular  Cu  pool  (fg/cell)  is  enhanced  as  fewer  ligands  are  available  to  keep  Cu 
from  partitioning  to  the  cell. 
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nM  Added  +  Ambient  Copper 


Figure  44.  Comparison  of  Thallasiosira  growth  in  0.4  pm  filtered  and  1 
kDa  permeate  samples  from  Norfolk  Site  A  (Hampton  Roads). 


Figure  44  above  shows  a  comparison  of  Thallasiosira  growth  in  0.4  pm  filtered 
and  1  kD  ultrafiltration  permeate  samples  from  Norfolk  Site  A  (entrance  to 
Hampton  Roads)  in  July  2001.  The  figure  clearly  illustrates  the  protective 
effect  of  colloidal  ligands.  Thallasiosira  grown  in  the  0.4  pm  filtered  sample 
(containing  the  complete  suite  of  complexing  ligands)  shows  no  reduction  in 
growth  until  a  threshold  of  at  least  300  nM  total  Cu  is  reached.  However,  in 
the  1  kD  permeate  (where  a  large  fraction  of  Cu  ligands  have  been  removed)  the 
toxic  effect  of  the  Cu  is  evident  at  a  much  lower  total  Cu  concentration  of  about 
120  nM. 
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Figure  45  below  presents  a  comparison  of  Cu  uptake  (fg  per  cell)  between  the 
0.4  pm  filtrate  and  lkD  ultra-filtration  permeate  from  Norfolk  Site  C  (Elizabeth 
River  end  member,  July  2001).  Uptake  patterns  are  remarkably  linear,  and  the 
sharp  contrast  in  slopes  is  consistent  with  relative  ligand  levels  ( 1  kD  permeate 
ligand  levels  much  lower,  allowing  greater  uptake  of  Cu  into  the  algal  cells). 


Figure  45.  Comparison  of  Cu  uptake  from  filtered  and  ultra-filtered 
fractions  at  NRF-C  (South  Elizabeth  River). 


nM  Added  +  Ambient  Copper 


A  San  Diego  Bay  transect  from  Bay  Mouth  to  Far  South  Bay  is  presented  in 
Figure  46.  Shown  are  filterable  (<0.4  pm)  Cu  concentrations,  levels  of  strong 
Cu-binding  ligand,  and  concentrations  of  Cu  associated  with  our  bioassay 
organism  (Thallasiosira) .  Despite  higher  Cu  levels  in  South  Bay,  uptake  of  Cu 
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into  Thallasiosira  is  low,  likely  reflecting  the  much  greater  abundance  of  strong 
Cu-binding  ligands  in  South  Bay.  However,  at  the  mid-Bay  Coronado  Bridge 
site,  ligand  levels  are  substantially  reduced  but  Cu  levels  remain  relatively 
high,  and  as  a  result  Cu  uptake  by  Thallasiosira  is  significantly  enhanced. 
Similarly  for  the  Bay  Mouth  site;  despite  quite  low  Cu  levels,  the  concentration 
of  bioavailable  Cu  is  actually  greater  than  in  South  Bay  as  strong  ligand  levels 
are  lower  than  filterable  Cu  levels. 


Figure  46.  Uptake  of  Cu  by  Thallasiosira  as  a  function  of  filterable  Cu 
concentrations  and  strong  Cu-binding  ligand  levels  along  a  San  Diego  Bay 
transect. 
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Table  4  presents  a  summary  of  bioavailability  measures  at  selected  field  sites 
as  a  means  to  show  a  few  key  findings  of  this  aspect  of  the  project. 
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Measurements  of  bioavailable  copper  at  the  contrasting  study  sites  indicate 
that  (1)  the  fraction  of  filterable  Cu  that  is  bioavailable  is  greatest  in  the  San 
Diego  system,  and  (2)  Cu  bioavailability  is  greatly  enhanced  by  ultrafiltration, 
where  up  to  75%  of  ambient  Cu  is  made  accessible  after  ligand  level  reduction. 


Table  4.  Bioavailable  Cu  fractions  at  selected  field  sites. 
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Figure  47  on  the  following  page  illustrates  how  these  data  can  be  placed  into  a 
modeling  framework  with  significant  predictive  power.  The  working  hypothesis 
is  that  the  levels  of  strong  Cu-binding  ligand,  LI,  (in  contrast  to  weaker  L2  and 
L3  species)  regulate  metal  availability  to  biological  systems  (e.g.  algal  cells). 
Therefore,  when  levels  of  total  filterable  Cu  are  less  than  LI,  uptake  into  cells 
should  be  prevented  or  at  least  hindered.  However,  when  Cu  levels  approach 
or  exceed  LI,  uptake  of  metals  by  organisms  should  be  relatively  unhindered 
and  cellular  metal  levels  should  rapidly  increase.  Indeed  this  is  what  we  are 
observing  in  the  three  study  systems  (refer  to  Figure  47).  If  we  plot  Cu  in  the 
probe  organism  (in  this  case  Thalassiosira,  fg  per  cell)  versus  excess  Cu  (i.e. 
[Cu]  >  [LI]),  we  observe  a  relationship  that  supports  the  working  hypothesis. 
Cellular  levels  of  Cu  are  extremely  low  until  [LI]  is  saturated,  at  which  point 
Cu  uptake  begins  to  increase  dramatically.  Remember  that  the  data  points  in 
this  plot  encompass  three  very  different  environments  (our  three  impacted 
harbors),  so  the  predictive  power  is  quite  robust.  Note  that  if  instead  of  the 
strong  ligand  [Li],  the  concentration  of  weak  ligand  [L2  or  L3]  is  substituted  in 
the  above  equation;  a  much  poorer  relationship  is  found,  therefore  reinforcing 
the  concept  of  strong  ligand  control. 
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Figure  47.  Copper  uptake  into  cells  can  be  accurately  predicted  from 
levels  of  strong  ligand. 
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XI.  High  Precision  Cu  Stable  Isotope  Ratio  Measurements 

A.  Goal  and  Summary 

Our  objective  was  to  determine  whether  sources  of  Copper  (Cu)  (and  in  later 
research,  Zinc  (Zn))  in  DOD-impacted  harbors  can  be  identified  by  unique 
stable  isotopic  signatures  of  these  transition  metals.  Reconstruction  and 
apportionment  of  historic  metal  inputs  may  also  be  possible  from  sediment 
cores  if  appropriate  coring  locations  can  be  identified.  Our  proof-of-principal 
evaluation  was  conducted  within  the  context  of  our  much  larger  study  of  Cu 
and  Zn  speciation  and  bioavailability.  When  we  initially  proposed  this 
hypothesis  and  effort,  we  felt  that  there  was  a  reasonable  chance  of  success 
and  in  general  the  potentially  large  benefits  (both  in  application  deliverables  for 
the  DoD  and  fundamental  geochemical  processes)  outweighed  the  risks. 
However,  it  was  certainly  true  that  this  was  cutting-edge  research  given  that, 
except  for  the  work  of  Marechal  et  al.  (1999)  and  Gale  et  al.  (1999),  that  few 
investigators  had  even  the  analytical  capabilities  to  quantify  the  isotope 
fractionation,  let  alone  the  environmental  application.  Our  group  possessed 
the  analytical  tools,  as  well  as  extensive  field  study  experience,  and  therefore 
was  in  a  unique  position  to  test  the  hypothesis.  We  initially  estimated  that  it 
would  require  18-24  months  to  reach  preliminary  conclusions,  and  as  it  turned 
out,  that  was  close,  but  somewhat  underestimated.  In  fact  it  took  us  roughly 
26  months  to  overcome  some  quite  recalcitrant  analytical  hurdles  and 
demonstrate  the  technical  feasibility  of  measuring  extremely  small  stable 
isotopic  ratio  variations  of  copper  in  some  very  nasty  matrices.  This  effort  was 
an  absolute  prerequisite  in  the  research  path  for  determining  the  viability  of 
using  copper  isotopic  fractionation  in  the  environment  to  trace  sources  of  the 
metal.  As  briefly  outlined  in  the  following  sections,  the  first  12  months  effort 
involved  ground-breaking  work  on  Cu  isotope  systematics  on  the  multi¬ 
collector,  magnetic  sector,  inductively-coupled  plasma  mass  spectrometer  (MC- 
MS-ICP-MS;  Micromass  Isoprobe).  The  bottom-line  conclusion  from  these 
studies  was  very  encouraging,  and  that  from  an  instrumental  viewpoint,  the 
applied  goal  was  instrumentally  technically  feasible.  We  then  moved  on  to 
developing  methods  for  recovery  and  cleanup  of  Cu  from  complex  matrices 
(marine  and  estuarine  waters  and  sediments)  This  work  has  been  particularly 
challenging  given  that  both  isolation/ extraction/ concentration  and  cleanup 
steps  must  be  essentially  quantitative  (i.e.  they  must  maintain  the  exact 
isotopic  signature  of  the  source  material).  It  was  also  determined  that  a  very 
high  degree  of  cleanup  was  required  because  the  MC-MS-ICP-MS  Cu  isotope 
ratio  precision  and  accuracy  was  significantly  influenced  by  the  presence  of 
other  cations/ anions  in  the  samples.  Though  this  might  sound  like  a  fairly 
straight  forward  exercise  in  good  analytical  chemistry,  the  added  constraint  of 
maintaining  (and  demonstrating)  endogenous  isotope  ratios,  as  it  turned-out, 
was  very  demanding  of  time  and  resources.  In  the  end,  however,  we  believe 
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that  we  have  overcome  the  technical/ analytical  hurdles  and  feel  that  a  pilot- 
scale  field  effort  is  certainly  justified.  We  worked  towards  that  goal  with  field 
samples  collected  in  parallel  with  the  speciation/ bioavailability  samples,  and 
initial  data  from  DoD-impacted  samples  showed  Cu  isotopic  signatures 
significantly  different  than  that  of  reference  standards  and  seawater  (Figure 
49).  That  anti-fouling  agent-impacted  samples  should  have  a  different  isotope 
signature  is  not  implausible  given  the  wide  range  of  Cu-isotope  fractionation 
now  recognized  to  exist  in  nature,  and  the  extreme  precision  of  our  methods. 
In  fact  it’s  very  gratifying  that  after  sticking  our  neck  out  several  years  ago, 
that  recent  reports  in  the  literature  confirm  the  early  work  of  Marechal  et  al. 
(1999)  and  suggest  that  Cu-isotope  fractionation  is  likely  to  prove  to  be  a 
powerful  new  tool  for  probing  geochemical  processes  and  source  attribution 
(Bermin  et  al.  (2006);  Chapman  et  al.  (2006);  Markl  et  al.  (2006);  Marechal  and 
Albarede.  (2002);  Zhu  et.  al,  (2002);  Coplen,  et.  al,  (2001);  Zhu,  et.  al,  (2000). 

Though  certainly  a  high  visibility  component  of  the  overall  project,  the  actual 
planned  fiscal  commitment  to  the  Cu-isotope  component  of  the  project  was 
quite  small  in  comparison  to  that  required  for  the  speciation  and  bioavailability 
components.  In  fact,  because  of  the  analytical  challenges  (now  resolved),  the 
actual  effort  devoted  to  the  Cu-isotope  work  went  well  beyond  what  was 
planned,  and  as  a  result  funding  and  time  for  a  comprehensive  evaluation  of 
Cu  isotope  ratios  in  all  the  field  samples  collected  was  inadequate.  However, 
the  hard  development  was  completed  successfully,  and  the  viability  of  the 
hypothesis  established.  Additional  work  will  be  required  to  field  validate  the 
extensive  isotope  effort  to  date.  We  believe  that  with  just  a  few  minor 
exceptions  that  field  estuarine  water  sites  are  adequately  sampled  to  date,  and 
that  these  should  provide  a  great  data  set  representative  of  in-harbor  mixed 
isotope  conditions.  In  a  couple  areas  we  feel  the  data  set  could  be 
strengthened: 

1 .  Copper  Ore  Samples 

We  have  analyzed  common  major  Cu-containing  reagents,  and  although 
Cu  isotope  fractionation  is  apparent  (which  bodes  well  for  actual  ore 
samples),  the  provenance  of  the  reagents  is  unknown.  In  collaboration 
with  the  Geology  and  Geophysics  Departments  we  have  obtained  a  small 
selection  of  common  ores,  and  these  remain  to  be  analyzed.  However,  it 
is  also  critical  that  we  identify  and  sample  the  current  and  historic 
feedstock  ores  for  common  Naval  and  Recreational  anti-fouling  coatings. 
We  believe  that  this  should  be  possible  using  contacts  within  the  Navy, 
Coating  Manufacturers,  and  Economic  Geologist  communities. 

2.  Anti-fouling  Paint  Samples 

In  addition  to  obtaining  a  more  complete  collection  of  “off-the-shelf’  anti¬ 
fouling  paints  representative  of  both  Navy  and  Recreational  boat  current 
and  historic  application,  we  would  approach  paint  manufacturers  for 
semi-processed  materials. 

Page  68  of  79 


68 


3.  Other  Navy  Cu-Sources 

Though  likely  a  relatively  minor  contribution,  Naval  ship  discharges  that 
leach  copper  from  internal  plumbing  and  equipment,  have  not  been 
adequately  sampled.  It  would  be  important  to  assess  whether  these 
sources  have  any  unusual  isotopic  signature  that  might  be  reflected  in  a 
disproportionate  impact  upon  isotopic  distribution  in  water  columns  or 
sediments. 

4.  Sediment  Cores 

Funding  and  time  did  not  permit  the  identification  and  sampling  of 
certain  selected  sediment  columns.  As  discussed  in  the  original 
proposal,  this  type  of  historic  record  would  be  invaluable  for 
demonstrating  and  comparing  the  impact  of  different  sources  of  Cu  to 
the  “impacted”  harbors. 


B.  Background  and  Theory 

Recent  advances  in  inorganic  mass  spectrometry  instrumentation  likely  will 
enable  a  novel  and  extremely  powerful  approach  to  Cu  and  Zn  source 
apportionment  -  the  use  of  stable  isotopic  signatures  of  these  metals  to 
deconvolute  sources.  While  source  apportionment  has  been  widely  applied 
with  Pb  (Rosman  et  al.,  1993;  Chiaradia  et  al.,  1997;  Sturges  and  Barrie,  1987) 
and  Sr  (Dia  et  al.,  1992;  Yoshihiro  et  al.,  1999;  Goldstein  and  Jacobsen,  1987) 
stable  isotopic  systems,  the  extreme  method  precision  required  for  Cu  and  Zn 
application  has  until  recently  prevented  a  systematic  evaluation  of  this 
potentially  valuable  method. 

The  nominal  stable  isotopic  distribution  of  Cu  and  Zn  in  the  terrestrial 
environment  is  tabulated  below: 


Mass  (Isotope) 

Copper (%) 

Zinc  (%) 

63Cu  62.9296  (63) 

69.09 

64Zn  63.9291  (64) 

48.89 

65Cu  64.9278  (65) 

30.91 

66Zn  65.9260  (66) 

27.81 

67Zn  66.9271  (67) 

4.11 

68Zn  67.9248  (68) 

18.56 

70Zn  69.9253  (70) 

0.62 
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Recent  studies  by  Chapman  et  al.  (2006);  Markl  et  al.  (2006);  Zhu  et  al.  (2002), 
Zhu  et  al.  (2000),  Marechal  et  al.  (1999)  and  Gale  et  al.  (1999)  and  earlier  work 
by  Shields  et  al.  (1965)  and  Walker  et  al.  (1958),  have  shown  that  stable  Cu 
isotope  ratios  in  copper  ores  vary  significantly  in  nature,  and  that  the  potential 
therefore  exists  for  source  discrimination  based  upon  isotopic  signature. 
Isotopic  variation  of  over  9  per  mil  between  various  ore  bodies  has  been 
reported  in  these  studies.  Copper  is  relatively  mobile  in  geologic  processes, 
and  through  incomplete  redox  and  exchange  reactions,  isotopic  fractionation 
may  occur. 

In  the  earliest  published  paper  addressing  natural  relative  isotopic  variation  in 
Cu,  Walker  et  al.  (1958),  using  a  Nier  type  mass  spectrometer,  measured  a  total 
spread  of  -1  to  +8  per  mil  in  a  number  of  minerals.  Measurement  error  was 
estimated  to  be  less  than  0.1%  using  an  NBS  CuO  primary  standard  and  a  Cul 
working  standard.  Samples  (n=12)  of  chalcocite  and  native  copper  from  the 
White  Pine  mine  in  Ontonagon  Michigan  (USA),  and  samples  (n=ll)  of 
volborthite,  malachite,  and  azurite  from  the  Cougar  mine  in  Montrose  Colorado 
(USA)  were  studied  along  with  16  other  samples  (including  brochantite, 
pitchblende,  marine  sediments,  niccolite).  The  sediment  sample  with  the 
highest  positive  enrichment  contained  a  large  amount  of  organic  matter  and 
was  collected  in  an  area  with  significant  anaerobic  bacterial  activity.  The  later 
finding  suggests  that  biologically-mediated  redox  cycling  may  enhance  isotopic 
fractionation. 

The  most  comprehensive  study  to  date  of  natural  variations  in  the  isotopic 
composition  of  Cu,  was  conducted  by  Shields  et  al.  (1965).  Over  100  samples, 
including  a  large  selection  of  primary  and  secondary  Cu-containing  minerals 
(chalcopyrite,  bornite,  chalcocite,  oxides,  carbonates,  native  Cu,  tetrahedrite, 
etc.)  were  studied.  These  investigators  employed  a  much  improved  mass 
spectrometer  compared  with  that  used  by  Walker  et  al.(1958),  and  were  able  to 
obtain  stable  ratios,  essentially  independent  of  time  as  long  as  the  ion  current 
was  not  decaying.  Shields  et  al.  (1965)  referenced  their  ratios  to  a  standard 
different  from  that  used  by  Walker  et  al.  (1958),  and  therefore  the  absolute 
fractionation  reported  (-9.0  to  +3.3  per  mil)  is  not  directly  comparable  with 
Walker’s  data  (conversion  is  easy  however),  but  the  total  range  in  fractionation 
observed  (12.3  per  mil)  is  remarkably  consistent.  The  authors  concluded  that 
variations  in  Cu  isotope  fractionation  appear  to  be  related  to  conditions  under 
which  a  given  mineral  species  forms,  with  a  secondary  influence  of  location. 

Marechal  et.  al.  (1999)  employed  modern  plasma-source  mass  spectrometry  to 
the  study  of  Cu  and  Zn  isotopic  fractionation  in  silicates,  ores,  and  biological 
materials.  These  investigators  referenced  their  ratio  data  to  the  NIST  Cu 
standard  (SRM  976),  a  convention  that  we  followed  in  our  work.  This  study  is 


Page  70  of  79 


especially  relevant  from  the  analytical  standpoint,  in  that  our  analysis  (both 
instrumental  and  clean-up)  protocols  incorporate  and  build  upon  the  detailed 
work  of  these  investigators.  Ratio  precision  of  20-40  ppm  (0.04  per  mil)  was 
reported,  similar  to  that  revealed  on  our  instrument  in  preliminary  validation 
studies.  Remarkably,  these  investigators  also  reported  a  range  in  Cu  isotope 
fractionation  of  Cu  ores  of  about  9  per  mil.  The  smaller  number  of  ores 
examined  (n=9)  prevented  any  definitive  conclusion  as  to  the  relationship 
between  mineral  type  and  isotopic  signature. 

The  work  of  Gale  et  al.  (1999)  is  important  in  at  least  three  major  areas:  (a) 
they  corroborated  earlier  studies  by  demonstrating  that  stable  isotope 
fractionation  of  Cu  in  natural  Cu-containing  minerals  range  between  -1.63  to 
+7.71  per  mil  (b)  an  exploratory  study  of  the  effect  of  smelting  and  fire  refining 
on  Cu  isotopes,  suggested  that  isotopic  fractionation  does  not  occur  in  these 
processes,  therefore  preserving  the  original  isotopic  signature  of  the  ore.  (c) 
they  compared  traditional  TIMS  techniques  with  state-of-the-art  MC-MS-ICP- 
MS  methods  (our  method  of  choice)  and  found  comparable  natural 
fractionation. 

The  important  initial  question  to  for  us  to  address  is  whether  the  background 
geochemical  isotopic  signatures  of  the  various  study  sites  are  significantly  and 
usefully  different  than  that  of  DoD  and  other  anthropogenic  copper  sources. 
For  source  apportionment  to  be  successful,  a  fractionation  of  >0.4  per  mil  must 
be  observed  between  major  end-members,  given  an  isotopic  ratio  precision  of 
20-40  ppm. 


C.  Analytical  Considerations 

Ultra-high-precision  isotope  ratio  measurements  of  Cu  and  Zn  require  the  use 
of  complex  and  expensive  mass-spectrometry  instrumentation.  Potentially  two 
classes  of  technologies  can  be  applied  to  the  analytical  challenge:  (A)  Thermal 
Ion  Mass  Spectrometry  (TIMS),  and  Multiple-Collector,  Magnetic  Sector, 
Inductively-Coupled  Plasma  Mass  Spectrometry  (MC-MS-ICP-MS).  The  former 
approach,  though  still  evolving  incrementally,  is  a  mature  technique,  having 
been  used  now  for  several  decades.  Its  primary  application  has  been  in 
systems  with  3  or  more  stable  isotopes,  where  the  ionization  potential  of  the 
element  is  relatively  low,  and  where,  due  to  radiogenic  parents,  large  natural 
variations  in  stable  isotopic  abundance  are  found.  These  elements  include: 
Neodymium  (Nd),  Lead  (Pb),  and  Strontium  (Sr).  However,  application  of  TIMS 
to  Cu-isotope  systematics  is  problematic  in  that  the  Cu  ionization  potential  is 
relatively  high  leading  to  very  poor  efficiencies,  and  double-spike  techniques  to 
correct  for  instrumental  isotopic  fractionation  cannot  be  applied  since  only  two 
stable  isotopes  exist.  In  addition  the  TIMS  analysis  is  extremely  time- 
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consuming,  which  coupled  with  the  required  pre-analysis  chemical 
separations,  dramatically  limits  sample  throughput. 

Instrumentation  for  the  later  approach,  MC-MS-ICP-MS,  has  been 
commercially  available  for  only  about  8  years,  and  it’s  only  in  the  last  5  years 
that  as  the  scientific  community  recognized  the  potential  of  this  technique  and 
compelled  manufacturers  to  modify  designs  and  correct  deficiencies,  that  the 
early  promise  has  been  fulfilled.  The  great  advantages  of  this  technique 
emanate  from  the  argon  plasma-front  end.  The  inductively-coupled  plasma  is 
a  very  efficient  ionization  source,  and  therefore  elements  with  high  ionization 
potentials  can  be  effectively  measured.  Also  analysis  time  is  shortened  and 
sample  preparation  simplified  over  traditional  TIMS  approaches.  One 
drawback  to  this  method,  though  not  technical  in  nature,  is  the  high  cost  of 
the  instrumentation,  typically  over  $750,000.  The  multi-collector  design 
allows  for  simultaneous  collection  of  several  ion  beams,  therefore  overcoming 
the  precision  limitation  imposed  by  plasma  “flicker”  in  single  collector 
instruments.  Multi-collector  instruments  can  produce  remarkably  precise 
isotope  ratio  data,  in  certain  systems  better  than  0.0005%  (5  ppm,  .005  per 
mil) .  This  level  of  precision  is  more  than  adequate  to  expect  that  stable  isotopic 
fractionation  in  non-radiogenic  element  systems  can  be  accurately  quantified. 
Particularly  exciting  is  the  potential  for  1st  and  2nd  row  transition  metals 
(including  Cu),  where  multiple  redox  states  are  the  rule,  and  where  biological 
and  geochemical  redox  cycling  may  lead  to  significant  isotopic  fractionation. 


C.  Analytical  Progress 

The  nearly  1  million  dollar  cost  of  the  MC-MS-ICP-MS  system  purchased  by  the 
UW-Madison  in  early  2000  was  funded  through  non-DOD  sources.  Our 
instrument  is  of  a  unique  design,  where  the  electrostatic  focusing  sector  of 
traditional  double-focusing  designs  is  replaced  by  a  hexapole-collision  cell.  The 
hexapole  improves  transmission  over  competing  designs,  thermalizing  and 
focusing  the  ions,  while  the  integrated  collision  cell  destroys  polyatomic  species 
endemic  to  plasma-based  instrumentation,  and  which  would  otherwise 
interfere  with  analyte  mass  spectra.  The  later  factor  cannot  be  over 
emphasized,  spectral  interference  removal  is  critical  for  accurate  and 
intercomparable  data. 


The  first  year  of  effort  in  this  area  of  study  was  dedicated  to  demonstrating  the 
technical  feasibility  of  measuring  extremely  small  stable  isotopic  ratio 
variations  of  copper.  This  work  was  a  necessary  first  step  in  determining  the 
viability  of  using  copper  isotopic  fractionation  in  the  environment  to  trace 
sources  of  the  metal.  These  studies  involved  ground-breaking  work  on  Cu 
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isotope  systematics  on  the  multi-collector,  magnetic  sector,  inductively-coupled 
plasma  mass  spectrometer  (Micromass  Isoprobe).  As  mentioned  previously,  the 
bottom-line  conclusion  from  these  studies  was  very  encouraging,  and  that  from 
an  instrumental  viewpoint,  the  applied  goal  was  technically  feasible. 


Analysis  Number 


Figure  48.  Copper  Stable  Isotope  Ratio  Precision:  MS-MC-ICP-MS 

An  example  of  the  external  precision  obtained  during  one  of  our  standard 
validation  experiments  is  shown  above  in  Figure  48.  Remarkable  precision  is 
apparent,  with  variation  not  evident  until  the  4th  or  5th  decimal  place.  A  per 
mil  standard  deviation  (1  o)  of  0.022  is  calculated,  more  than  an  order-of- 
magnitude  better  than  would  be  required  to  resolve  the  unit  per  mil  variations 
measured  in  copper  ores  (see  above). 

An  important,  but  somewhat  cautionary  outcome  of  these  studies  was  the 
finding  that  Cu  isotope  ratio  precision  and  accuracy  was  significantly 
influenced  by  the  presence  of  other  cations  in  the  samples.  Therefore,  we 
chose  a  very  conservative  approach  in  confronting  the  next  phase  of  the 
validation,  i.e.  isolation  and  concentration  of  Cu  from  complex  natural 


Page  73  of  79 


matrices.  Considering  the  technical  complexities  of  the  overall  goal,  we  decided 
that  to  dispatch  any  concerns  (either  our  own,  or  from  external  reviewers) 
about  matrix  effects,  that  we  would  purify  samples  such  that  only  trace  levels 
of  interfering  ions  remained.  Though  this  might  sound  like  a  fairly  straight¬ 
forward  exercise  in  good  analytical  chemistry,  the  added  constraint  of 
maintaining  (and  demonstrating)  endogenous  isotope  ratios,  as  it  turned-out, 
was  more  of  a  challenge  than  expected.  Our  focus  for  much  of  the  remaining 
project  period  was  to  develop  methods  for  recovery  and  cleanup  of  Cu  from 
complex  matrices  without  isotopic  fractionation  (i.e.  maintaining  exact!!!!! 
63/65  ratio).  We  directed  a  large  portion  of  this  effort  towards  methods 
capable  of  isolating  metals  from  relatively  large  volumes  of  marine  and 
estuarine  water. 

a.  Extraction:  Isolation  of  Cu  and  Zn  from  Seawater  Matrix. 

To  allay  concerns  about  isotopic  fractionation,  the  method  must  be  essentially 
100%  efficient  in  isolating  Cu  and  Zn,  but  equally  important  one  must  also  be 
able  to  recovery  100%  of  the  Cu  and  Zn  from  the  “sorbant”.  Some  isotopic 
fractionation  is  essentially  unavoidable  so  quantitative  recovery  is  critical  and 
also  serves  as  a  good  predictor  of  method  effectiveness.  We  tried  several 
approaches  including  both  solid  and  liquid-liquid  extractions.  For  example: 
various  immobilized  ligands  (iminodiacetate  on  polymer  backbone  (Chelex), 
other  direct  solid-phase  extractant  resins  (Eichrom;  Diphonix  and 
Monophosphonic).  Chelex  displayed  good  efficiency  in  extraction  from 
seawater,  but  we  were  not  able  to  achieve  100%  removal  from  the  Chelex,  even 
with  an  array  of  common  eluants.  Eichrom  resins  showed  the  opposite 
problem  -  poor  and  irreproducible  extraction  efficiency  of  Cu  from  seawater, 
but  good  recovery  from  the  solid  phase.  All  these  methods  were  also  plagued  by 
poor  efficiency  in  pre-cleanup  (i.e.  interfering  ions  were  also  extracted  with 
target  Cu  and  Zn).  Somewhat  discouraged,  we  then  turned  to  a  method  first 
proposed  by  Wells  and  Bruland  (1998)  for  isolation  of  several  metals  from 
saline  waters.  We  adapted  and  modified  this  method  and  after  several  months 
of  experimentation  feel  confident  that  this  approach  will  meet  our  stringent 
requirements.  A  brief  summary  of  the  technique  follows. 

One  must  first  synthesize  the  metal-binding  ligand  (bis(2-hydroxyethyl) 
dithiocarbamate  (HEDC).  We  did  so  according  to  the  method  of  King  and  Fritz 
(1985)  using  ethanolamine,  sodium  hydroxide,  and  carbon  disulfide  dissolved 
in  methanol,  followed  by  careful  precipitation  of  the  crude  ligand  upon  addition 
of  2-propanol.  The  ligand  was  subsequently  purified  and  characterized.  The 
clean  ligand  is  added  to  seawater  samples  to  form  a  neutrally-charged,  polar 
(water  soluble)  metal  ligand  complex.  The  ligand  has  high  affinity  and 
specificity  for  Fe,  Co,  Ni,  Cu,  Zn,  and  Cd.  The  neutral  metal  complexes  are 
recovered  onto  a  polystyrene-based  C-18  resin  column.  The  metal  complexes 
are  eluted  from  the  column  with  an  acidic  methanol  solution.  The  eluates  are 
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evaporated  to  near  dryness  after  addition  of  nitric  acid  to  destroy  the  organic 
ligand,  and  the  residue  taken  up  in  7  N  HC1.  Near  quantitative  isolation  of  the 
target  metals  from  major  ion  matrix  is  achieved  in  this  method.  In  brief,  the 
major  steps  in  this  process  are: 

Quantitative  Extraction  of  Cu  from  Complex  Matrices 

1.  Synthesize  ligand  (HEDC) 

2.  Purify  and  characterize 

3.  Form  metal-ligand  complex 

4.  Recover  complex  on  C-18 

5.  Elute  metal  from  column 

6.  Destroy  residual  ligand 

7.  Take  up  metals  in  7  N  HC1 


b.  Clean-up:  Separation  of  Cu  from  Other  Trace  Elements: 

The  sample  extract  is  further  cleaned-up,  and  trace  elements  separated  by 
anion  exchange  chromatography  on  AG  MP-1  resins.  Copper  is  eluted  from  the 
resin  with  23  mL  of  7  N  HC1.  Iron  and  Zn  can  be  subsequently  recovered  with 
2  N  HC1  and  0.5  N  HNO3  respectively.  This  clean-up  and  separation  method 
was  adapted  from  work  published  by  Marechal  et  al.  1999.  The  HCL  eluates 
are  evaporated  to  dryness  and  residue  taken  up  in  1%  HNO3  for  analysis  by 
MC-MC-ICP-MS. 

Clean-up  and  Separation  of  Cu  from  other  Elements 

1.  Anion  exchange  (AG  MP-1)  chromatography  using  7  N  HC1 
2.  Collect  fractions,  dry  down,  take  up  in  dilute  nitric  acid 


E.  Samples  for  Stable  Isotope  Analysis 

a.  Bio-fouling  Panel  Studies 

We  collaborated  with  Dr.  Peter  Seligman  (SPAWAR/SYSCEN)  to  evaluate  the  Cu 
(and  other  metals)  leached  from  panels  treated  with  Cu-biofouling  paints.  Dr. 
Seligman  provided  us  a  series  of  seawater  leachate  samples  (BRA-640)  from 
which  we  worked  to  provide: 

i.  high  precision  Cu  isotope  ratio  measurements  (after  extraction  and 
removal  of  seawater  matrix). 

ii.  comprehensive  multi-element  measurements  on  panel  extracts.  These 
data  will  provide  a  very  valuable  data  set  on  minor  and  trace  element 
leaching  rates  from  the  panels.  In  addition  these  data  will  likely  support  our 
source  reconciliation  efforts. 
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b.  Harbor/ Estuary  Studies 

One  liter  and  0.5  liter  samples  for  stable  isotope  measurements  were  collected 
from  each  of  our  study  sites  on  every  sampling  campaign.  Since  this  sampling 
program  focused  on  the  filterable  phases,  all  isotope  samples  were  field  filtered 
through  0.4  pm  Meissner  polypropylene  filters.  We  have  demonstrated  that 
this  filter  behaves  ideally  with  respect  to  Cu  and  Zn,  neither  sorbing  “dissolved” 
metal  nor  introducing  metal  contamination.  This  behavior  is  a  prerequisite  for 
accurate  isotope  ratio  characterization  of  the  samples.  We  also  collected 
particulate  phases  on  a  1  pm  Teflon  membrane  for  exploratory  measurements 
of  stable  isotope  ratios  in  suspended  particles.  Water  samples  were  stabilized 
with  ultra-high  purity  nitric  acid  to  ensure  high  metal  recovery.  The  blank 
contributed  by  the  acid  is  negligible. 

Impacted-Harbor  Cu-Isotope  Sample  Summary  (number  of  samples  collected) 

Cape  Fear 

Dissolved  (4) 

Particulate  (4) 

Colloidal  (6) 

Norfolk  Harbor 

Dissolved  (12) 

Particulate  (12) 

Colloidal  (5) 

San  Diego  Bay 

Dissolved  (28) 

Particulate  (28) 

Colloidal  (12) 

In  addition  to  the  Harbor/ Estuary  samples,  we  obtained  the  following  “source” 
samples. 


Cu  Stable  Isotope  Source  Sample  Summary 

Navy  ship  hull  panel  leachates  (BRA-640) 

A  commercial  Cu-based  antifouling  paint 

Six  Cu  ores  representing  industrial  feed  stocks 

Reference  seawater 

Vegetation  from  Cu  mine-impacted  area 

Several  rivers  draining  contrasting  watersheds 

Figure  49  on  the  following  page  presents  some  preliminary  Cu-isotope 
fractionation  (signature)  data.  Though  the  data  set  is  still  small,  and  we’re  not 
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yet  capable  of  asserting  that  our  hypothesis  is  definitive,  it  is  apparent  that  a 
large  range  in  isotopic  fractionation  is  exists  (more  than  an  order  of  magnitude 
greater  than  our  analytical  precision)  and  that  impacted-harbor  samples  are 
differentiated  from  selected  source  samples. 


White  Pine  Mine 
Copper  Salts 
Copper  Metal 
Seawater 
NBS  976 
BRA-640 
Norfolk  FD 
Norfolk  FC 
San  Diego  NF04 
San  Diego  F03 

-0.5  0.0  0.5  1.0  1.5  2.0  2.5  3.0 

del65Cu  (vsHPSCu) 


o 

o 

o 

Copper  Isotope 
Fractionation 
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Figure  49.  Measured  copper  stable  isotope  ratios  in  selected  reference 
and  field  samples. 
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